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Introduction 


1.1  Problem  Scenario 

The  expanding  technology  base  exploiting  millimeter  and  near-milli¬ 
meter  electromagnetic  waves  has  greatly  enhanced  the  ability  of  re¬ 
mote  sensors  to  detect  and  identify  complex  targets  of  interest.  The 
enhancement  is  especially  significant  when  the  target  is  obscured  by 
fog,  clouds,  or  smoke  which  can  render  infrared  and  optically  based 
sensors  useless.  As  a  result,  a  large  number  of  sensors  operating  at 
millimeter  and  near-millimeter  wavelengths  have  been  produced  for 
such  applications  as  identification  of  targets,  terminal  guidance,  and 
fuzing.  Accordingly,  the  need  for  efficient  iind  informative  simulation 
techniques  for  these  high-frequency  sensors  has  also  risen.  Consider 
the  scenario  depicted  in  figure  1. 

A  source/receiver  (S/R)  moving  along  some  trajectory  illuminates 
a  target  of  interest.  It  is  desired  to  estimate  the  return  from  the  target 
as  the  S/R  moves  along  the  trajectory.  Notice  that  whether  figure  1 
describes  a  tsuget  detection/identification  problem  or  the  terminal 
phase  of  a  guided  munition  is  mostly  determined  by  the  trajectory 
being  considered.  A  backscatter  simulation  technique  which  places 
few  or  no  restrictions  on  the  paths  of  the  trajectories  to  be  simulated 
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Figure  1:  Target  encounter  scenario. 
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could  therefore  find  use  in  all  phases  of  seeker  munitions  studies.  In 
addition,  of  course,  the  relative  motion  between  the  S/R  and  the 
target  could  be  due  strictly  to  the  motion  of  the  target.  Thus  the 
scenario  could  sdso  include  the  return  from  passing  targets.  It  will 
be  assumed  that  the  return  should  be  calculated  about  a  thousand 
times  per  meter  of  trajectory.  The  reason  for  mentioning  this  detadl 
here  is  to  introduce  early  the  impact  of  the  trajectory  on  any  practical 
modeling  method.  In  particular,  this  study  attempts  to  include  the 
trajectory  into  the  structure  of  the  problem,  as  opposed  to  simply 
considering  the  trajectory  as  a  sequence  of  points  at  which  to  evaluate 
the  return. 

This  report  begins  with  a  brief  review  of  a  modeling  method  devel¬ 
oped  at  the  Harry  Diamond  Laboratories  for  the  computer  simulation 
of  fuze/air-target  encounters.  The  simulation  technique  uses  a  com¬ 
posite  of  simple  analytic  structures  (cones,  spheres,  flat  plates,  etc) 
to  model  the  surface  of  a  target  of  interest.  A  geometric  optics  as¬ 
sumption  is  then  used  to  estimate  the  expected  amplitude  and  phase 
of  the  return  as  the  S/R  moves  along  the  trajectory  of  interest.  The 
current  study  makes  use  of  bicubic  B-spline  patches  to  describe  ar¬ 
bitrary  complex  targets  in  greater  detail.  Using  the  geometric  optics 
approach  for  estimating  the  target  return,  the  backscatter  problem 
becomes  one  of  finding  adl  the  specular  points  of  such  an  arbitrary 
surface  with  respect  to  the  S/R  location. 

The  main  portion  of  this  report  describes  a  method  for  exploiting 
the  continuity  of  the  resulting  surface  model  to  reduce  the  problem 
of  finding  all  the  specular  points  associated  with  each  new  trajectory 
position  to  that  of  tracking  the  motion  of  existing  points.  In  particu- 
leir,  it  is  shown  that  the  locations  of  the  annihilations  and  creations  of 
speciilar  paths  may  be  predicted  for  an  entire  trajectory,  eliminating 
the  need  to  seeirch  the  whole  surface  for  specular  points  as  the  target 
moves  in  relation  to  the  S/R.  It  is  further  shown  that  the  same  tech¬ 
nique  can  be  extended  to  include  so-called  multiple-bounce  return,  in 
which  it  is  assumed  that  the  signal  is  reflected  about  the  surface  two  or 
more  times  before  returning  to  the  receiver.  The  general  flow  of  sim¬ 
ple  algorithms  for  the  implementation  of  the  new  strategy  is  described 
and  the  results  of  example  runs  displayed.  The  examples  point  out 
potential  robustness  problems  with  the  simplistic  algorithms,  prob¬ 
lems  which  are  not  inherent  in  the  modeling  method  itself.  The  final 
section  of  the  report  outlines  potential  solutions  to  these  problems  as 
well  as  extensions  to  the  model.  The  potential  solutions  and  exten- 


sions  are  presented  as  part  of  the  suggested  areas  of  further  reseairch 
and  include 

•  the  potential  for  further  exploitation  of  the  properties  of  bicu¬ 
bic  B-spline  surfaces  to  increase  the  robustness  of  the  specular 
tracking  algorithms, 

•  extensions  to  include  surface  discontinuities, 

•  extensions  to  include  multiple  bodies  and  shadowing,  and 

•  experimental  validation  and  comparison  with  other  methods. 

A  Useful  Backscatter  Modeling  Technique 

Methods  for  developing  radar  backscatter  models  of  air  targets  suit¬ 
able  for  computer  simulation  of  radar- fuze/ mr-target  encounters  axe 
described  by  Dammann  [1].  The  models  are  developed  in  three  steps: 

First,  the  target  is  approximated  by  an  ensemble  of 
simple  geometric  shapes  such  as  ogives,  cylinders,  flat 
plates,  and  ellipsoids.  Second,  the  specular  point  ( the  point 
where  the  surface  is  perpendicular  to  the  incident  rays)  is 
located  on  each  shape.  All  the  return  is  assumed  to  come 
from  these  specular  points,  and  the  strength  of  the  return 
from  each  point  is  assumed  to  depend  on  the  local  curva¬ 
tures  of  the  surface  at  that  point.  Third,  the  returns  from 
all  specular  points  are  summed  vectorially  to  yield  the  total 
target  return.  [1] 

This  modeling  method  has  proved  to  be  very  useful  both  for  pre¬ 
dicting  the  backscatter  characteristics  of  interest  in  fuze  design  and 
for  the  simulation  of  fuze  processors.  In  particular,  because  the  model 
provides  actual  Doppler  voltage  waveforms,  the  computer-generated 
waveforms  can  be  fed  directly  into  actual  fuze  processors  to  predict 
their  response. 

The  target  modeling  technique  described  above  has  been  particu¬ 
larly  successful  in  the  air-target  encounter  scenario  partially  because 
of  the  relative  simplicity  of  the  shapes  of  most  aircraft.  In  many  cases, 
between  10  and  20  primitive  analytic  shapes  are  all  that  are  necessary 
to  capture  the  essential  shape  of  a  particular  aircraft.  As  the  targets 
of  interest  become  more  complex  or  the  desire  to  match  their  surface 
more  accurately  increases,  the  use  of  simple  analytical  shapes  to  de¬ 
scribe  the  tauget  surface  often  becomes  impractical.  Tensor-product 


piecewise  polynomial  surfaces  are  a  simple  means  of  describing  com¬ 
plex  surfaces,  and  as  such,  present  a  practical  way  to  extend  the 
geometric  optics  approach  to  include  targets  of  increased  complexity. 

1.3  Surface  Modeling  with  B-splines 

In  this  paper  it  is  assumed  that  a  three-dimensional  vector  function 
R(«,t;)  exists  which  describes  the  target  surface  of  interest;  that  is, 
given  a  particular  choice  of  u  and  v,  R(u,  v)  represents  the  i-,  y-,  and 
^-coordinates  of  a  point  on  the  target  surface.  It  is  further  assumed 
that  over  the  entire  domain  of  R  (call  this  U  xV)  the  partial  derivative 
vectors 

du'dvi  ^  ’ 

exist  and  are  continuous  for  i  +  j  <  2.  B-splines  can  be  used  to  con¬ 
vert  a  patchwork  of  data  points  into  just  such  a  function.  For  the 
major  residts  of  this  paper  it  is  the  second-order  continuity  of  the 
surface  which  is  exploited  to  reduce  the  backscatter  simulation  time, 
and  an  understanding  of  B-splines  is  really  not  necessary.  However, 
the  following  brief  discussion  may  be  helpful.  (The  reader  interested 
in  the  details  of  B-splines  in  general  is  encouraged  to  study  de  Boor’s 
book  A  Practical  Guide  to  Splines  [2]  or,  for  a  quicker  look,  see  Gordon 
and  Riesenfeld  [3],  Cox  [4],  or  Amos  [5].) 

The  mathematical  term  “spline”  originates  from  the  mechanical 
spline  used  by  draftsmen  to  trace  out  smooth  curves.  When  thin 
strips  of  wood  or  metal  are  pinned  at  strategic  locations  on  the  draft¬ 
ing  board,  the  material  takes  on  a  shape  to  minimize  its  internal  stress 
and  gives  a  generally  pleasing  interpolant  to  the  data.  If  the  drafts¬ 
man  is  not  pleased  with  the  resulting  trace,  the  number  of  pinning 
points  and  their  locations  may  be  altered  until  the  desired  shape  is 
formed.  Piecewise  polynomial  functions  are  the  mathematical  mate¬ 
rial  of  choice  for  fitting  smooth  curves  through  or  near  a  given  set  of 
data.  The  term  “piecewise”  refers  to  breaddng  up  the  entire  curve  into 
sections  of  low-order  polynomials  (most  often  cubic).  The  sections  are 
joined  at  “knots”  so  that  appropriate  continuity  is  mzdntained.  The 
term  “B-spline”  refers  to  a  special  set  of  piecewise  polynomial  func¬ 
tions  which  form  a  basis  for  edl  piecewise  polynomials  having  a  given 
knot  set  and  continuity  constraints  at  the  knots.  As  an  example,  the 
cubic  B-spline  basis  functions  with  uniformly  spaced  knots  having  ze¬ 
roth,  first,  and  second  derivative  continuity  at  each  of  the  knots  may 
be  written  as 


8 


5i(z)  =  (l-z)V6, 

Biiz)  =  (3z=*  -  6z*  +  4)/6, 

Bsiz)  =  (-3z®  +  3z*  +  3z  +  l)/6, 

Bi(z)  =  2^/6. 

Note  that  a  linear  combination  of  the  functions  above  can  be  used  to 
construct  any  desired  cubic  polynomied.  Now  suppose  that  we  have 
a  rectangular  mesh  of  points  which  we  would  like  to  use  to  describe 
a  surface.  The  mesh  can  be  used  to  form  a  tensor  product  surface 
using  the  same  basis  functions.  Consider  the  behavior  of  the  following 
function: 

4  4 

.=1  i=i 

where  Ai+i„j+i  are  three-dimensional  vectors  of  finite  length  and  u  and 
V  are  real  numbers  ranging  from  0  to  1.  v)  is  simply  one  way  of 

writing  a  bicubic  B-spline  surface  defined  over  a  uniform  rectangular 
knot  set  for  the  A:,/th  patch  of  the  quilt  of  data  points.  Clearly,  for 
any  fixed  vedues  of  k  and  /,  and  for  all  u  and  v  6  (0, 1),  the  various 
partial  derivatives  of  equation  (1)  are  continous  functions  of  u  and  v 
for  all  positive  integers  i  and  j.  It  is  also  easy  to  verify  that  all  the 
zeroth,  first,  and  second  parti£d  derivatives  with  respect  to  u  and  v 
are  continuous  across  neighboring  values  of  k  and  1.  For  example, 

d'^Rk,i{u,v)  __  ^^dBi{u)^  ^  ^dBj[v) 

dudv  du  dv 

=  1/4 1-(1  -  u)*[— (1  —  vY Ak+\,i+i  +  —  4t;)Afc+i,(+2 

-t-(— 3u*  -f  2t;  4-  l)Afc+i,i+3  -f  u*Afc+i_j+4] 

-|-(3u*  —  4«)[— (1  —  v)*Afc+2,j+i  +  (3v*  —  4v)Afc+2./+2 
-f-(  — 3v*  +  2v  +  l)Afc+2,t+3  +  V*Aj|t+2,/+4] 

-t-(— 3u*  -f  2u  -b  1)[— (1  —  v)*Afc+3,i+i  -t-  (3v*  —  4v)Afc+3,n.2 
-b(— 3v*  -t-  2t;  -b  l)Afc+3,j+3  -b  v*A|i,+3,j+4] 

+(ti*)[-(l  —  u)*Afc+4.j+i  +  (3v*  -  4v)Afc+4,{+2 
-b(  — 3u*  -b  2w  -b  1)A|H.4,J+3  -b  V*Afc+4,j+4]|  , 
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from  wliich  it  is  an  easy  check  that 


1/4  {—(1  —  n)*[Afc+i,j+i  +  Afc+i^j+a] 

+(3'u*  —  4u)[Afc+2,j+i  +  Afc+2,J+3] 

+  (— 3u*  +  2u  -f  1)[A*+3,/+i  +  A*4.3,/+3] 
+(u*)[Afc+4,J+l  +  Afc+4,/+3]} 

dudv 

While  many  surface  modeling  methods  offer  smooth  surfaces,  the  B- 
spline  surface  in  generzd  has  three  major  properties  which  make  it 
particularly  well  suited  for  surface  modeling.  The  first  is  that  the  sur¬ 
face  passes  near  (in  a  sense  that  will  not  be  elaborated  here)  the  points 
A,j.  This  allows  the  modeler  to  choose  a  mesh  of  points  which  gen¬ 
erally  describes  the  desired  surface  (much  like  a  planar  facet  model) 
and  then  observe  the  resulting  surface.  As  might  be  expected,  the 
more  points  that  are  used  to  define  the  surface,  the  closer  the  sur¬ 
face  will  pass  to  the  points.  The  second  major  property  is  the  local 
influence  of  the  points  A,j-;  that  is,  a  change  in  the  position  of  the 
point  Aij  only  changes  the  surface  near  the  point.  The  local  nature 
of  the  spline  allows  the  modeler  to  move  points  around  to  improve  the 
surface  in  a  particular  area  without  modifying  the  rest  of  the  model. 
The  third  major  property,  closely  related  to  the  first,  is  the  variation- 
diminishing  property  of  B-splines.  Loosely  speaking,  the  variation- 
diminishing  property  guarantees  that  any  plane  will  cut  through  a 
three-dimensionsil  B-spline  curve  no  more  times  than  it  does  a  linear 
interpolant  to  the  original  data  points. 

It  is  the  second-order  continuity  of  the  surface  which  is  exploited 
for  the  main  resiilts  of  this  paper.  Since  many  other  methods  exist  for 
representing  smooth  surfaces,  the  maiin  results  are  more  general  than 
just  for  B-spline  surfaces.  However,  mostly  because  of  the  previously 
stated  B-spline  properties,  the  surfaces  exemplified  in  this  report  are 
all  based  on  B-splines.  In  addition,  a  large  portion  of  the  suggested 
further  research  is  aimed  towards  exploiting  those  properties  specific 
to  B-spline  surfaces. 
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2.  Geometric  Optics  Approach  Using  B-splines 


2.1  Single-Bounce  Return 


The  geometry  of  the  specular  return  problem  from  a  single  patch 
of  an  arbitrary  B-spline  surface  is  shown  in  figure  2.  Rt(A)  is  the 
current  position  of  the  projectile  along  a  lincM  trajectory,  R,(u,i>) 
describes  the  target  surface  as  a  function  of  the  two  parameters  u  and 
17,  and  G(u,  v,  A)  is  the  difference  between  the  two  vectors  Rt(  A)  and 
R  ,{u,v).  It  is  assumed  throughout  the  analysis  that  the  projectile  ha^ 
an  unobstructed  view  of  the  surface  being  considered  (the  problem  of 
shadowing  is  ignored  and  addressed  only  in  the  section  on  suggested 
further  research).  A  necessary  and  sufficient  condition  for  a  point  on 
the  surface  to  be  a  specular  point  relative  to  the  position  Rt(  A)  is  that 
the  I2  norm  of  G(u,  17,  A)  be  either  a  local  maximum  or  minimum  with 
respect  to  the  two  target  surface  parameters  u  and  v.  Finding  all  the 
specular  points  of  a  given  surface  (for  a  given  trajectory  position)  is 
therefore  the  same  as  finding  all  u  and  v  that  satisfy  the  two  nonlinear 
equations 


and 


Ou 


dv 


=  0 

(2) 

=  0. 

(3) 

If  R(u,v)  is  given  by  a  tensor  product  of  cubic  B-splines  on  a  uni¬ 
form  grid,  both  Fi(u,v,A)  and  F2{u,v,A)  may  be  written  explicitly 
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in  terms  of  u  and  v.  However,  the  u  and  v  which  solve  equations  (2) 
and  (3)  (c^lll  them  u*  and  v*)  cannot  be  found  explicitly,  since  Fi(u,v) 
is  a  bivariate  polynomial  of  degree  five  in  u  and  six  in  v,  while  v) 
is  a  bivariate  polynomial  of  degree  six  in  u  and  five  in  v.  In  generzd 
then,  solving  equations  (2)  and  (3)  for  u*  and  v*  requires  a  numer¬ 
ical  lechnique.  Simple  application  of  Newton’s  method  for  nonlinear 
equations  will  find  solutions  to  (2)  and  (3)  provided  the  search  is  be¬ 
gun  “close  enough”  to  The  question  of  how  close  is  close 

enough  is  a  complex  one  which  ultimately  depends  on  the  variation 
of  the  surface  being  considered.  An  explanation  of  Newton’s  method 
and  algorithms  which  incorporate  the  method  may  be  found  in  most 
optimization  texts  including  Luenberger  [6].  The  essentiad  step  is  to 
solve  the  following  linear  system  of  equations: 

J(u,v)S  =  F(u,  v),  (4) 

where  F(u,u)  is  the  2x1  vector  [Fi(ti,u),F2(u,u)]^,  J(u,u)  is  the 
2x2  Jacobian  of  F(«,u)  with  respect  to  u  and  v,  and  S  is  the  2x1 
vector  for  the  new  step  in  u  and  v.  A  derivation  of  the  Newton  step 
for  the  location  of  a  specular  point  on  a  B-spline  surface  may  be  found 
in  appendix  A.  It  is  easy  to  see  that  for  the  bicubic  B-spline  surface 
being  considered,  the  elements  of  J(u,u)  are  continuous  functions  of 
u  and  u;  that  is,  all  second-order  derivatives  of  G(u,v)  exist  and  are 
continuous  functions  of  u  and  v. 

Using  Newton’s  method  alone,  it  would  be  possible  to  thoroughly 
search  each  B-spline  patch  for  specular  points,  keeping  track  of  their 
positions  and  local  principal  radii  of  curvatures  for  use  in  the  geomet¬ 
ric  optics  estimation  of  the  return  from  the  composite  model.  The 
procedure  would  have  to  be  repeated  for  each  new  position  filong  the 
trajectory.  To  get  a  feel  for  the  time  which  may  be  required  for  a  sin¬ 
gle  simulation,  suppose  a  target  were  being  modeled  by  B-splines  over 
an  80  by  50  knot  grid  (4000  patches,  give  or  take  a  few  end  conditions) 
and  that  it  were  possible  to  thoroughly  search  100  B-spline  patches 
per  second.  The  specular  return  from  the  target  could  be  calculated 
in  approximately  40  s  for  one  point  along  the  trajectory  of  interest. 
Now  suppose  further  that  we  are  interested  in  a  10-m  trajectory  to 
be  sampled  at  1-mm  increments  (not  an  unlikely  requirement  when 
working  with  millimeter  waves).  The  simulation  of  the  signature  for 
the  single  trajectory  would  take  approximately  400,000  s,  or  a  little 
over  four  and  one  half  days.  Admittedly,  other  techniques  could  prob¬ 
ably  be  used  to  cut  down  on  the  time  required  to  search  the  composite 
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model  for  specular  points.  (Methods  based  on  recursive  subdivision  of 
the  tcirget  surface  using,  for  example,  the  Oslo  algorithm  [7]  should  be 
considered.)  However,  any  method  which  searches  all  4000  patches  for 
each  position  of  the  trajectory  is  going  to  be  very  computer  intensive. 


2.2  Specular  Paths 


Whether  gradient  techniques  such  as  Newton’s  method  are  fast  or  not, 
they  can  give  important  insight  into  the  motion  of  specular  points.  Let 
the  current  trajectory  position  (Rt)  of  figure  2  vary  linearly  between 
the  two  endpoints  Rti  and  Rt2  as 


Rt(A)  —  ARt2  +  (1  —  A)Rtj. 


(5) 


For  the  surface  considered  here,  the  distance  function  G(u,v,A)  is 
a  continuously  differentiable  function  of  «,  v,  and  A.  Suppose  the 
coordinates  of  a  specular  point  are  known  for  a  paurticular  value  of 
A  and  we  wish  to  observe  the  motion  of  the  specular  point  as  A 
changes.  The  following  argument  is  a  trivial  extension  of  that  given 
by  Longuet- Higgins  for  a  time- varying  anaJytic  surface  [8].  The  con¬ 
ditions  for  a  specular  point  are  given  by  equations  (2)  and  (3).  Taking 
the  differential  of  equations  (2)  and  (3)  with  respect  to  u,  v,  and  A 
yields 


du^ 

yilG(..,.,,A)f 

dudv 


dudv 


dud^ 


^  y||G(u  ,  A)r^^  ^ 

dv^  dvdA 


which  may  be  rewritten  in  matrix  form  as 


J(u,v,A) 


du/dA 

dv/dA 

a»IIG(u,v.A)ii» 

auask 

0»||G(«.v,A)||» 

dvdA 


(6) 


where  J(u,v,A)  is  the  same  2x2  Jacobian  used  in  the  Newton’s 
method  search  (with  its  dependence  on  A  explicitly  noted).  Since 
the  elements  of  J(w,v,  A)  are  continuous  functions  of  u,  v,  and  A,  if 
J(u,u,  A)  is  nonsingular,  du/dA  and  dvfd^  will  both  be  finite,  which 
implies  that  the  changes  in  u  and  v  can  be  kept  as  small  as  desired 
by  choosing  the  change  in  A  small  enough.  Longuet- Higgins  [8]  refers 
to  the  vanishing  of  the  determinant  of  J(w,u,  A)  as  a  twinkle.  The 
physical  significance  of  this  result  is  that  as  the  source/receiver  moves 


across  a  second-order  continuous  surface,  specular  points  cannot  sud¬ 
denly  appear  or  disappear  unless  A)  is  locally  singular.  The 

observation  that  specular  points  move  in  continuous  paths  broken 
only  when  A)  is  singular  leads  directly  to  the  following  conclu¬ 

sion.  If,  for  any  given  trajectory,  it  were  possible  to  determine  all  the 
points  {u,v,A)  for  which  J(u,v,  A)  is  singular,  it  would  no  longer  be 
necessary  to  search  the  entire  surface  for  specular  points  at  different 
positions  along  the  trajectory.  It  would  only  be  necessary  to  find  all 
the  specular  points  corresponding  to  one  trajectory  position  (Rti  for 
example)  and  then  track  their  motion,  picking  up  or  losing  specular 
paths  only  at  twinkles.  It  is  shown  in  section  2.4  that  the  annihila¬ 
tion  or  creation  of  specular  paths  must  also  occur  in  pairs;  that  is,  a 
pair  of  specular  paths  must  run  into  each  other  to  be  annihilated,  or 
originate  from  a  single  point  on  the  surface  to  be  created. 

To  get  a  feel  for  the  potential  computational  savings  to  be  gained 
by  first  finding  the  twinkles  associated  with  a  particular  trajectory, 
consider  the  previous  example.  Suppose  that  on  the  average  there 
were  10  specular  points  on  our  surface  of  4000  B-spline  patches  and 
that  it  took  about  the  same  amount  of  time  to  update  the  location 
of  5  specular  points  eis  it  did  to  search  an  entire  patch  for  specular 
points.  Once  all  the  twinkles  have  been  found  for  the  trajectory  and 
the  problem  is  reduced  to  a  specular  tracking  problem,  it  will  take 
on  the  average  1/2000  times  as  long  to  find  each  new  set  of  specular 
points;  thus,  the  computational  time  will  be  reduced  from  a  little  over 
4-1/2  days  to  3-1/3  minutes,  plus  whatever  time  it  takes  to  find  the 
twinkles.  Recalling  that  the  search  for  twinkles  is  only  conducted 
once  for  a  given  trajectory,  we  could  reduce  the  computationeil  time 
required  to  simulate  the  backscatter  from  complex  targets  by  two  or 
more  orders  of  magnitude. 

The  following  example  both  clarifies  how  twinkles  may  be  used  to 
trace  the  motion  of  specular  points  along  a  surface  for  a  given  trajec¬ 
tory,  and  motivates  the  derivations  that  follow.  Figures  3  through  8 
show  the  evolution  of  a  specular  plot  for  a  simple  tzu-get  surface  and 
an  associated  trajectory  of  interest.  To  keep  the  example  relatively 
uncluttered,  it  is  assumed  that  in  addition  to  having  the  continuity 
properties  discussed  before,  the  surface  is  such  that  no  specular  points 
lie  on  its  edge  for  any  location  along  the  trajectory.  The  first  step  is 
to  find  all  the  specular  points  associated  with  the  beginning  of  the 
trajectory;  in  this  example,  only  one  such  point  exists  on  the  surface, 
as  shown  by  the  X  in  figure  3.  The  next  step  is  to  find  all  the  twinkles 
associated  with  the  chosen  trajectory  and  partition  the  trajectory  ac- 


14 


Figure  6:  After  birth  of  two  specular  paths 
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!  ;  cordingly.  In  addition  it  will  be  helpful  to  know  whether  each  of  the 

Hi  twinkles  represents  a  birth  or  a  death  of  a  pair  of  specular  paths.  It 

\  has  already  been  shown  that  at  a  twinkle  the  Jacobian  which  might 

ii  ordinarily  be  used  to  determine  the  location  of  the  specular  points 

r  becomes  singular,  making  it  difficult  to  locate  specular  paths  near  a 

twinkle.  Fortunately,  a  solution  is  available  to  the  problem  of  predict- 

■  ing  the  motion  of  specular  points  near  a  twinkle,  as  well  as  predicting 

'  whether  a  twinkle  represents  the  birth  or  death  of  a  pair  of  specu¬ 

lar  paths,  eis  shown  in  section  2.4.  Figure  4  shows  the  partitioning 
of  the  trajectory,  labeling  of  the  twinkles  according  to  whether  they 
represent  births  or  deaths,  and  prediction  of  specular  paths  near  the 

p  twinkles.  For  all  trajectory  locations  between  the  origin  of  the  tra- 

;  jectory  and  the  first  twinkle,  it  is  only  necessary  to  determine  the 

I  :  motion  of  the  original  specular  point  (since  no  new  specular  points 

L  will  occur  emy where  on  the  surface  except  at  a  twinkle).  Figure  5 

[  shows  the  motion  of  the  original  specular  point  in  the  region  preced- 

■  ing  the  first  twinkle.  When  a  twinkle  representing  the  birth  of  two 

I  specular  paths  is  reached,  the  two  new  paths  must  be  included  in  the 

I  set  of  specular  points  to  be  tracked,  at  least  until  the  next  twinkle. 

I  Figure  6  shows  the  situation  after  the  first  twinkle  and  shortly  before 

L  the  second.  At  that  time  one  of  the  specular  points  which  originated 

P  at  the  birthing  twinkle  was  heading  for  an  annihilation  with  the  orig- 

j  inal  speculw  point.  Figure  7  shows  the  specular  plot  at  a  time  after 

I  the  two  specular  paths  collided  at  the  remaining  twinkle.  For  the  re- 

I  mainder  of  the  trajectory,  it  is  only  necessary  to  track  the  motion  of 

L  the  one  remaining  specular  path,  as  shown  in  figure  8.  Once  again,  it 

*  should  be  pointed  out  that  the  14  by  14  grid  of  spline  patches  which 

formed  the  surface  in  this  example  was  only  searched  once  for  the  ini- 

I  tial  specular  points  and  then  once  again  for  the  twinkles.  After  that, 

1  specular  points  were  simply  tracked. 

I  2.3  Finding  Twinkles 

I  As  shown  in  the  previous  section,  a  necessary  condition  for  the  cre- 

J  ation  or  smnihilation  of  a  specular  path  is  the  vanishing  of  the  determi- 

M  nant  of  J(u,v,  A).  Denote  the  determinant  of  J(u,i;,  A)  as  F3(ti,v,  A). 

*  Then  F3(u,v,  A)  is  given  by 

1 

1  a>G(«,«,A)8>G(»,v,A)  /'9>G(»,»,A)\' 

F3K»,a)  =  — - - (  -  j  ■ 
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In  addition  to  F3{u,v,  A)  vanishing,  the  conditions  of  a  specular  point 
must  still  be  satisfied  in  order  for  the  point  (u,v,A)  to  be  a  twinkle. 
Given  a  target  surface  and  a  particular  trajectory  of  interest,  the 
twinkles  for  the  trajectory  may  be  found  by  solving  the  nonlinear 
equations 

Fi(ti,v,  A)  =  0,  (8) 

F2(u,v,A)  =  0,  (9) 

F3(u,v,A)  =  0,  (10) 

for  the  three  unknowns  u,  v,  and  A.  Once  again,  simple  gradient 
techniques  may  be  used  to  search  for  the  solutions  of  (8)-(10).  The 
actual  Newton  steps  for  finding  the  twinkle  coordinates  u,  v,  and  A 
on  a  uniform  B-spline  mesh  are  derived  in  appendix  B. 

Birthing  Specular  Paths 

By  themselves,  the  conditions  for  a  twinkle  do  not  tell  whether  a 
particular  twinkle  represents  a  birth  or  death  of  a  pair  of  specular 
points  with  respect  to  a  chosen  positive  trajectory  direction.  Once 
again,  a  straightforward  extension  of  the  analysis  by  Longuet-Higgins 
[8]  gives  a  method  for  describing  the  motion  of  speculzur  paths  near  a 
twinkle,  including  whether  the  twinkle  represents  the  birth  or  death 
of  a  pair  of  paths.  Consider  the  following  substitutions  in  Longuet- 
Higgins  analysis  for  the  motion  of  specular  points  near  a  twinkle: 


f(x,y,0 


=  A, 

=  |iG(«,u,A)|^ 

||G(u,u,A)||^ 

dti'dv^ldA'^ 


u=v=a=o 


where  the  coordinate  system  is  chosen  such  that  u  =  u  =  A  =  0at 
the  twinkle  and 

flOOO  =  fllOO  =  <Jk)10  =  Olio  =  0200  =  0.  (11) 

It  can  be  easily  seen  from  Longuet-Higgins  analysis  that  near  a  twinkle 
the  u,  V  coordinates  of  a  specular  point  Me  given  by 

r— 2aioi  A] 

u  =  ±  -  , 

O300 

O210O101  —  O300O011 

V  =  - ^ - ,  (12) 

O300O020 
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Equations  (12)  show  that  if  aioi/usoo  is  positive,  then  two  solutions 
exist  when  A  is  less  than  zero  and  no  solutions  exist  when  A  is  greater 
than  zero;  that  is,  a  pair  of  speculau^  paths  is  annihilated.  Similarly,  a 
pair  of  speculu  paths  is  created  when  a\oi/ 0,300  is  negative.  It  remzdns 
to  determine  a  transformation  of  coordinates  for  which  equations  (12) 
hold.  Simply  choosing  the  origin  of  the  new  coordinate  system  to  be 
the  location  of  the  twinkle  ensures  that  oooo  is  equal  to  zero.  It  will 
be  useful  here  and  for  analysis  to  come  to  consider  the  surface  formed 
by  letting  the  function  (iG(tt,r,  A)||*  be  the  w  coordinate  in  an  or¬ 
thogonal  u,  V,  w  coordinate  frame,  as  shown  in  figure  9.  For  lack  of 
a  better  term,  this  surface  is  referred  to  as  the  “distance  surface.” 
The  condition  for  a  twinkle  may  then  be  interpreted  physically  as  the 
vanishing  of  the  Gaussian  curvature  of  the  distance  surface.  Or  alter¬ 
nately  stated,  at  a  twinkle,  one  of  the  two  principal  radii  of  curvature 
of  the  distance  surface  is  equal  to  zero.  The  curvature  of  any  smooth 
surface,  at  a  given  point,  in  the  direction  6u,  6v  may  be  written  as 
(see,  among  others.  Faux  and  Pratt  [9]) 

Kn  =  (^u)*a200  +  (^tt)(^w)aiio  +  (^v)*a020.  (13) 

Now  suppose  the  coordinates  are  rotated  so  that  the  u-coordinate  is 
aligned  with  the  principal  radii  of  curvature  having  zero  value  (such 
a  direction  must  exist  at  a  twinkle).  Then  with  6v  =  0,  equation  (13) 
yields 

0  =  (^«)*a2oo> 


Figure  9:  Distance  surface. 
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which  implies  that  0200  =  0.  It  can  be  seen  that  the  twinkle  condition 
implies 

02000020  ~  Ojio  =  0. 

Therefore,  in  the  rotated  coordinate  system,  ono  must  also  be  equal 
to  zero  and  equations  (11)  are  satisfied.  Define  the  new  coordinates 
u',  and  A'  by 

u  =  o'cos^  —  u' sin^ -I- Uf„, 

V  =  u  sin6  +  v' cosO  +  vtwy 
A  =  A'  +  «t„, 

where  u^u,,  Vtwy  and  are  the  original  coordinates  of  the  twinkle. 
Then  if  $  is  the  angle  which  rotates  the  original  u-axis  into  the  di¬ 
rection  of  zero  curvature,  the  signs  of  oioi  and  O300,  in  the  new  coor¬ 
dinates,  will  tell  whether  the  twinkle  represents  a  birth  or  a  death. 
In  addition,  the  motion  of  the  specular  points  near  the  twinkle  (in 
the  new  coordinates)  will  be  given  by  (12).  The  derivation  of  the 
equations  for  0300,  Uioi,  <*210)  non,  and  0020  in  terms  of  the  original 
coordinates  and  6  is  given  in  appendix  C. 

2.5  Rotation  of  Coordinates 


In  order  to  carry  out  the  procedure  of  the  preceding  section,  it  is  nec¬ 
essary  to  determine  the  angle  B  which  rotates  the  original  u-direction 
into  the  direction  of  zero  (principal)  radius  of  curvature.  Define 
U(t)  :=  (u(t),v(t))  to  be  a  curve  on  the  surface  R(u(<),T;(t)).  It 
can  be  shown  (see  for  example.  Faux  and  Pratt  [9],  p  112)  that  for 
dU  fdt  to  be  in  the  direction  of  a  principal  radius  of  curvature  it  must 
satisfy  the  following  matrix  equation: 


where 


[D  -  /c„E](iU/df  =  0 


[D]  = 


^^rCw,v) 


g^r(u,«) 
“ ■  Sugv 
0^r(u,p)  g*r(u,v) 

“  ■  dytSv  “  ■ 


> 


gr(u,v)  ^  dr{u,v)  ^r(u.v)  ,  Or(u.v) 

du  du  du  dv 

&t(u^v)  ^  Bt(u,v)  dt(u,v)  ^  ' 

du  dv  dv  dv 

and  n  is  the  local  unit  normal  to  the  surface.  The  conditions  that 
dU/dt  point  in  the  direction  of  a  principal  radius  of  curvature  and 
that  the  curvature  be  zero  require  that 


[D](nJ/d<  =  0. 


(14) 
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Ordinarily  matrix  equation  (14)  would  have  only  the  trivial  solution 
dU /dt  =  0.  However,  for  the  case  being  considered  here,  the  elements 
of  [D]  are  the  same  eis  the  elements  of  the  Jacobian  introduced  earlier, 
which  has  been  shown  to  vanish  at  a  twinkle.  Therefore,  the  two 
equations  implied  by  (14)  are  not  independent,  and  we  are  free  to 
choose  either  as  the  necessary  condition  for  dU/dt  to  point  in  the 
direction  of  zero  curvature.  The  resulting  relation  is 

u/v  =  -dij/dii 
or 

B  =  t&n~^{-dn/dn) 

where  6  is  shown  in  figure  10  and  dij  is  the  t,jth  element  of  [D]. 
Note  that  it  does  not  matter  in  which  of  the  two  directions  of  zero 
curvature  we  choose  to  rotate  the  u-eods,  since  specular  points  will  be 
annihilated  or  created  in  both  directions. 

Once  a  twinkle  is  found,  it  is  a  straightforward  matter  then  to 
C2dculate  the  angle  6  which  rotates  the  original  u-axis  into  the  direc¬ 
tion  of  zero  (principal)  radius  of  curvature,  using  the  first  row  of  the 
Jacobian.  The  equations  derived  in  appendix  B  can  then  be  used  to 
find  the  values  of  aaoo,  Uioi,  a2io,  ooii»  atnd  0020  in  the  new  coordi¬ 
nate  system.  Comparison  of  the  signs  of  O300  and  ajoi  tells  whether 
the  associated  twinkle  represents  a  birth  or  a  death  of  two  specular 
paths.  In  addition,  equations  (12)  can  be  used  to  approximate  the 
motion  of  the  specular  points  near  the  twinkle.  The  approximation 
can  be  importsmt  in  the  actual  implementation  of  any  algorithms, 
since  the  vanishing  of  the  Jacobian  at  the  twinkle  makes  it  difficult  to 
use  ordinary  gradient  techniques  to  search  for  the  new  speculeir  paths 
indicated  by  a  birth  twinkle. 


Figure  10:  Rotation  angle  for  zero  curvature. 
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3.  Extension  to  Multiple  Bounces 

3.1  nth  Order  Specular  Points 

Often  multiple-bounce  return,  as  depicted  in  figure  11,  produces  a  sig¬ 
nificant  contribution  to  the  overall  target  backscatter.  Recall  that  the 
geometric  optics  approach  assumes  that  the  energy  ememating  from 
the  source  strikes  the  target  and  is  then  reflected  in  the  same  way  as 
would  a  plane  wave  incident  on  a  plane  tangent  to  the  surface,  or  put 
even  more  intuitively,  as  a  ball  bouncing  off  a  billiard  cushion.  For 
the  single-bounce  return  already  considered,  this  implies  that  energy 
will  only  return  from  locations  on  the  surface  for  which  the  local  sur¬ 
face  normal  points  towards  the  S/R.  However,  for  complex  surfaces, 
that  energy  which  is  reflected  by  the  surface  in  other  directions  may 
strike  another  surface  location  so  as  to  be  directed  back  to  the  re¬ 
ceiver.  When  the  magnitude  of  the  energy  returned  after  multiple 
bounces  is  calculated,  both  the  increased  distance  traveled  and  the 
loss  at  each  reflection  serve  to  diminish  the  result.  For  this  reason, 
the  single-bounce  return  is  often  referred  to  as  the  primary  or  first- 
order  specular  return,  while  the  double-bounce  return  is  referred  to  as 
the  secondary  or  second-order  specular  return  (with  the  obvious  nota- 
tionsd  extension  to  higher  order  bounces).  The  actual  target  surface 
locations  at  which  the  reflections  are  considered  to  occur  are  denoted 
the  nth  order  specular  points.  For  many  targets  of  interest,  it  may 
be  assumed  that  the  multiple-bounce  return  is  negligible  compared  to 
the  single-bounce  return.  However,  if  the  target  surface  is  assumed  to 
be  highly  reflective,  the  local  curvatures  of  the  surfaces  at  the  specular 
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Figure  11:  Multiple-bounce  return. 
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points  small,  and  the  added  round  trip  distance  not  too  large,  then 
multiple  bounces  may  occur  which  have  return  energy  with  magnitude 
of  the  same  order  as  the  single  bounces. 


3.2  The  Bistatic  Bounce 


Fortunately,  it  can  be  eeisily  shown  that  the  conditions  for  am  nth 
order  specular  bounce  are  directly  related  to  the  total  distance  of  the 
line  segments  connecting  the  source  to  the  receiver  via  each  of  the  nth 
order  specxilar  points.  In  fact,  the  following  result  can  be  stated  as  am 
immediate  consequence  of  a  weak  form  of  Fermat’s  principle  of  optics. 
However,  a  simple  derivation  is  included  here  in  hopes  of  giving  the 
reader  an  intuitive  feel  for  the  stationary  conditions  of  the  multiple- 
bounce  problem.  Consider  the  case  of  two  rays  originating  from  a 
point  O  on  am  arbitrary  surface  (z,y, /(z,t/)),  ais  shown  in  figure  12. 
Choose  the  points  at  amd  B  at  {xi,y2,Zi)  to  be  points 

on  the  two  different  rays  so  that  the  distance  AO  equads  the  distance 
Ob.  Without  loss  of  generality,  we  may  choose  the  origin  of  the 
system  to  be  at  the  point  O.  Define  a  third  point  C  at  (azj,  ayj,  az2). 
By  varying  a,  notice  that  the  point  C  can  be  moved  to  any  position 
along  the  ray  OB.  The  distamce  AOC  is  given  by 

^AOC  =  \/(*-*i)*  +  (y-yi)*  +  (/(®,j/)-^i)* 

+  v/(x  -  ozj)*  +  (y  -  oy2)*  -I-  (/(*,y)  -  022)* 


with 


dx 


+ 


1/^X0 

^l^oc 


2(z  -  zi)  +  2(/(z,y)  - 
|2(z  -  azj)  +  2(/(z,y)  -  az2 )  ~ g~~ ~ 


Evaluating  the  partiad  derivative  at  z  =  y  =  /(x,y)  =  0  and  setting 
the  result  equad  to  zero  yields 


^I^AO 


xi  + 


dx 


+  aj  Doc 


Z2  +  Z2 


df{x,y) 


dx 


=  0. 


Noting  that  a/D-^  =  the  condition  that  the  distamce 

AOC  be  stationary  with  respect  to  z  implies 


^  -(gi  +  X2) 
dx  (zi  +  Z2) 


(15) 
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Figure  12:  Rays  emanating  from  a  surface. 


Repeating  for  the  partial  derivative  of  the  distance  with  respect 
yields 

^/(»>y)  ^  -(yi  +  Vi) 

dy  (z-i  -h  Z2) 


to  y 
(16) 


Note  that  equations  (15)  and  (16)  are  not  defined  when  Zi  =  —Z2- 
However  the  coordinate  system  could  always  be  rotated  so  that  zi 
is  not  equal  to  Z2,  allowing  us  to  assume  without  loss  of  generality 
that  Zi  does  not  equal  — zj.  Equations  (15)  and  (16)  imply  that  two 
vectors  lying  in  the  tangent  plane  are 


Ti  =  (zi  4-  22)1  -  (xi  +  X2)k 

and 

T2  =  {z\  4  Z2)j  ~  (yi  +  yi)^. 

And  from  the  above  assumption,  Tj  and  T2  form  a  basis  for  the 
tangent  plane  at  O. 

Now  consider  the  conditions  for  a  bistatic  bounce.  A  ray  leaving 
point  A  and  striking  the  surface  at  O  will  be  reflected  through  C  if 
and  only  if  it  pzisses  through  B.  Since  Djq  is  equal  to  D^,  any  ray 
starting  at  (xi,yi,zi)  and  striking  the  surface  at  O  will  pass  through 
the  point  (xj,y2,  Z2)  if  and  only  if  the  ray  Rm  (shown  in  fig.  13)  is  nor- 
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Figure  13:  Bistatic  bounce. 


mal  to  the  surface  tangent  plane  at  O.  The  ray  connects  the  origin 
with  the  point  {Xm,ymtZm)  which  lies  midway  between  {xi,yi,zi)  and 
If  Hm  is  orthogonal  to  both  Ti  and  T2  as  defined  above, 
then  Rrn  will  be  normal  to  the  surface  tangent  plane  at  O. 

Rm  =  1/2{(*1  +  Xi)!  +  (pj  +  1/2)3  +  (^1  +  ^2)k] 

implies 

<  Rm,Ti  >=  1/2[(X1  +  X2)(Zl  +  Z2)  -  (X]  +  *2)(zi  +  Z2)]  =  0, 

<  Rm,T2  >=  l/2[(yi  -I- y2)(zi  +  Z2)  -  (yi  -t-y2)(^i  +  ^2)]  =  0. 

Under  the  mild  assumptions  that  the  surface  is  continuous  and  dif¬ 
ferentiable  with  respect  to  both  x  and  y,  if  the  distance  >1  to  O  to  C 
is  stationary  with  respect  to  the  surface  parameters  at  O,  then  the 
conditions  for  a  bistatic  bounce  are  satisfied. 

3.3  The  nth  Order  Bounce 

Consider  the  implications  of  the  previous  bistatic  result  for  a  double 
bounce.  The  double  bounce  (in  which  a  signal  originates  at  1,  is 
bounced  from  2  to  3,  and  then  from  3  to  1,  as  shown  in  fig.  14)  may 
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Figure  14:  Double  bounce. 

be  viewed  as  the  simultaneous  result  of  two  bistatic  return  problems: 
i.e.,  from  1  to  2  to  3  and  from  2  to  3  to  1.  Both  of  the  bistatic 
cases  must  satisfy  their  respective  local  stationary  conditions,  and 
so  the  total  distance  12  4-  23  4-  31  must  be  stationary  with  respect 
to  the  four  surface  parameters  used  to  describe  the  two  reflecting 
surfaces.  The  argument  immediately  generalizes  to  any  number  of 
bounces.  Thus,  just  as  the  single-bounce  specular  points  on  a  surface 
may  be  found  by  searching  for  stationary  points  of  the  distance  from 
the  source/receiver  to  the  surface  (and  back  if  you  like),  multiple- 
bounce  paths  may  be  found  by  searching  for  stationary  points  of  the 
total  round-trip  distance.  Each  new  bounce  simply  introduces  two 
new  stationary  conditions  along  with  two  new  surface  parameters. 
The  same  Newton  technique  which  was  used  to  find  single-bounce 
specular  points  may  be  used  to  find  multiple-bounce  specular  points 
(recall  that  we  are  ignoring  for  now  the  problem  of  shadowing).  The 
stationary  conditions  may  be  stated  in  terms  of  the  distances  shown 
in  figure  11  as 

duj  dvj 

for  all  j  =  1,  ...,  n.  Noting  that 

ddj  _  ddj  ^  ^ 
duj  dvj 
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whenever  t  —  1  <  j  or  j  >  t  +  2,  we  get  the  2n  conditions 

d{di  4-  dj+i)  _  Q  _  ^(dj  +  dj+i) 
dui  ,  dvi 

Denote  by  <7,-.  Then,  taking  the  differential  of  the  previous  2n 

equations  with  respect  to  the  2n  surface  parameters  =  l,n) 

and  the  trajectory  parameter  A  yields 


/  d^G, 

ail*  auia«i 
a*Gi  a*G, 

duidvi  aiq 


\  ailiaCn  dvidvn 


5*G,  d*G,  \ 

* 

r  d*G, 

duidun  duidVn 

k 

duidA 

d*Gi  d*G, 

dundvi  dVidVn 

dA 

dvidA 

d*Gn  d*G» 

dvn 

d*Gj 

dUndVn  dvl  / 

■  dA  ■ 

dv„dA  . 

(17) 

The  same  argument  that  led  to  the  conclusion  that  specular  paths 
can  only  be  created  or  annihilated  at  twinkles  may  be  directly  ex¬ 
tended  to  include  multiple-bounce  twinkles;  i.e.,  the  vanishing  of  the 
determinant  of  the  2n  x  2n  Jacobian  of  equation  (17)  is  required  for 
a  discontinuous  motion  of  the  nth  order  specular  paths. 


3.4  Example  of  Savings  Using  Multiple^Bounce  Twinkles 

Consider  again  a  B-spline  surface  made  up  of  4000  patches:  now  con¬ 
sider  the  computational  time  required  to  compute  the  second-order 
(double-bounce)  return  from  the  target.  The  search  for  second-order 
specular  points  may  be  conducted  by  checking  all  pairs  of  patches 
for  the  stationary  conditions  described  earlier.  If  we  know  that  a  ray 
leaving  the  source  strikes  the  ith  patch,  then  strikes  the  jfth  patch, 
and  finally  returns  to  the  receiver,  there  is  no  need  to  search  for  a  ray 
leaving  the  source,  striking  the  jth  surface  and  then  the  ith  surface, 
and  then  returning  to  the  receiver.  Note  however,  that  it  may  be  pos¬ 
sible  for  a  multiple  bounce  to  occur  on  a  single  patch.  The  number  of 
patch  combinations  which  must  be  considered  for  each  location  along 
the  trajectory  is  therefore  simply  the  number  of  distinct  combinations 
of  the  patches  taken  n  at  a  time,  where  n  is  the  number  of  bounces. 
The  number  of  distinct  patch  combinations  which  must  be  considered, 
given  k  patches  and  n  bounces,  is  simply 
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For  k  »  n  the  number  of  patch  combinations  is  closely  approximated 
by 


Once  again,  suppose  that  there  are  on  the  average  10  double-bounce 
specular  paths  (20  double-bounce  specular  points)  on  the  surface  for 
any  single  location  along  a  given  trajectory,  and  that  it  takes  about 
the  same  amount  of  time  to  update  the  location  of  5  paths  as  it  does 
to  search  an  entire  pair  of  patches  for  an  unknown  number  of  double¬ 
bounce  paths.  Once  all  the  double-bounce  twinkles  have  been  found 
for  the  trajectory  and  the  problem  is  reduced  to  a  specular  tracking 
problem,  it  will  take  on  the  average  4  x  times  as  long  to  search 
the  entire  target  surface  for  double-bounce  specular  paths  as  it  will 
to  track  the  paths.  More  generally,  let  Tm  denote  the  average  time 
required  to  update  the  n-bounce  specular  paths  and  T,  be  the  average 
time  required  to  search  an  n-patch  combination  for  an  unknown  num¬ 
ber  of  specular  points.  Then,  once  the  twinkles  and  initial  n-bounce 
specular  paths  have  been  found,  it  will  take  Q  times  as  long  to  search 
the  entire  surface  for  specular  paths  as  it  will  to  update  the  paths, 
where  Q  is  given  by 


Of  course,  it  would  be  hoped  that  a  one-time  search  of  the  surface  for 
patch  combinations  that  could  not  be  involved  in  a  multiple  bounce 
for  any  trajectory  would  drastically  reduce  the  number  of  patches  to 
consider.  However,  on  a  complex  surface,  the  number  of  potentizd 
multiple-bounce  patch  combinations  is  still  likely  to  be  very  large  and 
it  is  clearly  worth  spending  some  time  searching  for  the  multiple- 
bounce  twinkles. 


4.  Algorithm  Development  and  Examples 

4.1  Introduction 

The  flow  of  programs  required  to  implement  the  twinkle  method  is 
outlined  in  this  section,  and  examples  of  the  results  of  some  simple 
Fortran  77  realizations  of  these  algorithms  are  given.  The  simple 
programs  make  the  (occasionally  incorrect)  assumption  that  at  most 
one  twinkle  and/or  initial  specular  point  will  be  present  on  a  single 
B-spline  patch.  In  addition  to  examples  where  the  simple  programs 
work  without  a  hitch,  some  examples  are  discussed  to  show  the  effects 
of  missing  either  a  twinkle  or  an  initial  specidar  point.  It  is  also  shown 
that  in  many  cases  the  results  of  tracking  the  specular  paths  can  lead 
to  discovering  twinkles  that  were  initially  missed. 

4.2  Program  Flow 

Figure  15  is  a  flow  chart  of  a  simplistic  program  which  finds  the  motion 
of  single-bounce  specular  points  edong  the  target  surface  as  the  S/R 
moves  along  the  trajectory.  Once  the  desired  B-spline  target  surface 
and  trajectory  have  been  selected,  the  search  for  twinkles  and  initial 
specular  points  begins.  The  method  used  to  search  for  twinkles  is 
simply  to  start  at  the  center  of  each  patch  and  the  center  of  the 
trajectory,  and  then  make  the  Newton  step  derived  in  appendix  B.  If 
the  steps  are  such  that  the  search  is  led  off  the  patch  or  trajectory, 
the  search  is  resumed  at  the  appropriate  edge.  If  the  search  continues 
to  leave  the  patch  or  trajectory,  it  is  assumed  that  the  patch  has  no 
twinkles  for  the  given  trajectory.  Such  an  approach  will  never  find 
more  than  one  twinkle  on  a  single  patch  and  is  not  even  guaranteed 
to  converge  to  a  twinkle  if  it  is  the  only  one  on  the  patch.  However, 
the  method  is  adequate  for  the  purposes  of  demonstrating  the  overall 
simidation  technique  and  is  fairly  robust.  For  example,  when  the 
patches  are  small  compared  to  the  overall  target,  most  trajectories 
are  such  that  no  single  patch  has  more  than  one  twinkle.  In  addition, 
when  a  single  twinkle  does  exist  on  a  patch,  the  method  nearly  always 
converges  to  the  twinkle.  The  terms  “nearly  always,”  “smaJl,”  and 
“most”  are  used  very  loosely  in  this  discussion  and  are  meant  only 
to  convey  the  message  that  even  when  such  a  simplistic  technique  for 
searching  for  the  twinkles  is  used,  the  method  often  works.  When 
a  twinkle  is  found,  the  analysis  of  sections  2.4  and  2.5  is  used  to 
calculate  the  coefficients  of  the  Taylor  series  expansion  of  the  motion 
of  the  specular  points  near  the  twinkle,  as  well  as  to  determine  whether 
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Figure  15:  Simple  flow  chart  of  algorithm  for  tracking  specular  points. 
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the  twinkle  represents  the  birth  or  death  of  a  pair  of  specular  tracks. 

The  seuch  for  the  specular  points  associated  with  the  beginning 
of  the  trajectory  is  conducted  in  a  manner  analogous  to  the  search 
for  twinkles.  For  each  patch,  the  initial  guess  for  the  location  of  a 
specular  point  relative  to  the  trajectory  origin  is  at  the  center  of  the 
patch.  If  the  search  exceeds  some  criterion  for  leaving  the  patch,  it  is 
assumed  that  no  specular  points  exist  on  the  patch  at  the  beginning 
of  the  trajectory.  If  a  specular  point  is  found,  it  is  assumed  to  be  the 
only  one  on  that  patch. 

Once  the  twinkles  and  initial  specular  points  are  found,  it  is  a  sim¬ 
ple  matter  to  sort  the  twinkles  according  to  their  trajectory  position. 
With  the  assurance  that  no  new  specular  paths  will  be  created  before 
the  next  twinkle  is  reached,  tracking  the  motion  of  the  specular  points 
becomes  very  fast.  The  current  program  simply  uses  the  last  known 
location  of  each  specular  point  as  the  starting  point  for  a  Newton 
search.  The  search  for  each  updated  position  generally  takes  only  a 
few  steps  to  converge.  When  the  next  twinkle  is  reached,  the  action  to 
be  taken  of  course  depends  on  whether  the  twinkle  was  found  to  be  a 
creation  or  annihilation.  For  an  annihilation,  a  check  of  the  distances 
from  the  twinkle  to  each  of  the  specular  points  quickly  shows  which 
two  specular  points  are  colliding.  The  two  points  which  are  closest  to 
the  twinkle  are  simply  dropped  from  the  list  of  specular  points  to  be 
tracked.  If  the  twinkle  represents  a  creation,  then  the  Taylor  coeffi¬ 
cients  for  the  motion  of  the  specular  points  near  the  twinkle  are  used 
to  find  the  locations  of  the  two  new  specular  paths  shortly  after  the 
twinkle. 

4.3  Another  Single-Bounce  Example 

Figures  16  and  17  show  another  more  complex  exsunple  of  the  use  of 
twinkles  for  tracking  single-bounce  specular  paths  on  a  B-spline  sur¬ 
face.  Figure  16  shows  the  B-spline  surface  control  mesh  and  desired 
trajectory,  as  well  as  the  locations  along  the  trajectory  where  specu¬ 
lar  path  discontinuities  are  expected  (based  on  a  search  for  twinkles). 
Each  twinkle  location  along  the  trajectory  is  labeled  B  or  D,  depend¬ 
ing  on  whether  the  twinkle  represents  the  birth  or  death  of  a  pair  of 
specular  paths,  and  lines  have  been  drawn  connecting  them  with  their 
associated  locations  on  the  target  surface.  In  addition,  the  predicted 
paths  of  the  specular  points  in  the  area  near  each  twinkle  are  shown. 
Figure  17  shows  the  results  after  tracking  the  specular  paths  for  1000 
locations  along  the  trajectory.  A  comparison  of  figures  16  and  17 
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shows  that  the  specular  paths  did  in  fact  remain  continuous  every¬ 
where  except  at  the  twinkles  and  moved  as  predicted  in  the  regions 
near  each  twinkle.  Because  completely  searching  the  entire  target  for 
specular  points  at  each  trajectory  location  was  unnecessary,  the  entire 
simulation  took  only  a  couple  of  minutes.  Even  if  the  global  search 
for  specular  points  could  be  reduced  to  10  s  per  trajectory  location, 
the  simulation  would  have  taken  over  2-1/2  hours  without  the  use  of 
twinkles. 

4.4  Missing  a  Twinkle 

Consider  the  effect  of  missing  a  twinkle  representing  the  death  of  two 
specular  paths.  With  no  modifications  to  the  flow  chart  shown  in 
figure  15,  the  program  would  continue  endlessly  searching  for  the  two 
speculsur  points  which  would  still  be  presumed  to  exist  after  reaching 
the  location  of  the  missed  twinkle.  However,  such  a  situation  (a  miss¬ 
ing  twinkle  representing  ^ln  annihilation)  can  easily  be  checked  for.  A 
quick  calculation  would  show  that  near  the  lost  path(s),  two  specular 
points  were  very  close  together.  The  program  could  conclude  that  a 
twinkle  representing  a  death  had  been  missed,  drop  the  two  specular 
points  from  the  list  to  be  tracked,  and  proceed.  A  missing  twinkle 
which  represents  the  death  of  two  specular  paths  is  therefore  easy  to 
handle  and  causes  no  great  concern  to  the  program.  For  the  case  of  a 
finely  incremented  trajectory,  the  tracking,  creation,  and  suinihilation 
portions  of  the  current  simple  program  are  each  efficient  and  reliable. 
Experience  has  shown  that  if  the  program  cannot  update  a  specular 
point  which  is  supposed  to  be  there,  it  invariably  turns  out  that  the 
specular  point  has  been  annihilated  at  a  twinkle  which  was  missed  in 
the  original  search.  In  fact,  twinkles  which  represent  annihilations  of 
specular  paths  could  be  considered  somewhat  redundamt  checks  of  the 
specular  search. 

Simple  questions  added  to  the  algorithms  can  provide  feedback 
which  corrects  for  most  missing  initial  specular  points  and  twinkles. 
The  basic  premise  is  that  errors  at  one  stage  of  the  specular  tracking 
will  show  up  as  inconsistencies  in  later  annihilations  or  creations.  For 
example,  if  a  twinkle  representing  an  annihilation  is  reached  which 
has  only  one  approaching  speculu  path,  the  Taylor  expansion  can  be 
used  to  find  the  missing  path  near  the  twinkle.  Once  the  missing  path 
has  been  found,  the  trajectory  may  be  backtracked  until  the  source 
of  the  missing  path  is  reached.  The  source  of  the  missing  path  may 
then  be  tested  to  see  if  it  resulted  from  a  missed  initial  specular  point 
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or  a  twinkle  representing  a  birth.  If  the  original  error  was  a  missed 
birthing  twinkle,  the  program  must  find  and  track  forward  the  other 
specular  path  which  must  also  have  been  missed.  Such  a  back-and- 
forth  tracking  scheme  can  be  effective  in  finding  and  correcting  for 
the  occasionally  missed  twinkle  or  initial  speculzur  point.  Note  for 
example,  that  in  the  case  shown  in  figures  16  and  17,  if  any  one  of  the 
initial  specular  points  or  twinkles  had  been  missed,  they  could  have 
been  found  and  corrected  for  with  the  above  scheme. 

If  it  were  not  for  the  fact  that  a  pair  of  specular  paths  may  be 
created  which  never  annihilate  with  any  other  paths  for  a  particvilar 
trajectory,  it  might  be  possible  to  skip  the  search  for  twinkles  alto¬ 
gether.  However,  those  pairs  of  paths  that  were  created  at  a  twinkle 
and  either  annihilated  with  each  other  later  or  never  amnihilated  at 
all  wotild  not  be  found  under  such  a  scheme.  It  should  also  be  pointed 
out  that  the  use  of  feedback  to  search  for  missing  twinkles  or  initial 
specular  paths  requires  a  history  of  edl  the  specular  paths  to  be  kept, 
so  that  the  formerly  lost  paths  may  be  added  in  appropriately.  Note 
that  if  no  initial  specular  points  or  twinkles  are  missed,  then  the  simu¬ 
lated  backscatter  may  be  calculated  at  each  step  and  there  is  no  need 
to  record  the  actual  paths.  Fortrsin  implementations  of  the  feedback 
strategy  have  been  successfully  developed.  However,  while  feedback 
may  be  used  to  correct  for  most  missing  twinkles  and  initial  speculsir 
points,  it  is  still  clearly  to  om  advantage  to  pursue  methods  which 
reduce  the  chances  of  either. 

4.5  Double-Bounce  Example 

Figures  18  and  19  show  an  example  of  the  use  of  double-bounce  twin¬ 
kles  for  tracking  double-bounce  speculjir  points  on  a  simple  B-spline 
surface.  Figure  18  shows  a  control  mesh  for  a  simple  crescent-shaped 
ribbon  which  has  been  tilted  slightly  so  it  may  be  viewed.  A  short 
trajectory  is  shown  along  with  the  only  double-bounce  twinkle  associ¬ 
ated  with  that  trajectory  and  sm^ace.  In  addition,  a  triangle  has  been 
drawn  to  show  the  double-bounce  path  associated  with  the  twinkle. 
Figure  19  shows  the  results  of  searching  for  double-bounce  specular 
points  at  40  locations  along  the  trajectory  (starting  at  the  lower  left 
of  the  trajectory).  As  expected,  the  number  of  double-bounce  specu¬ 
lar  points  associated  with  each  trajectory  location  before  the  twinkle 
did  not  chsuige  (there  were  none).  At  the  twinkle,  two  sets  of  specu¬ 
lar  paths  were  created  (pairs  of  specular  bounces  associated  with  the 
same  double-bounce  path  are  shown  connected  by  a  line).  The  two 
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Proposed  Additional  Research 


5.1  Introduction 

The  use  of  twinkles  to  reduce  the  computationzil  time  required  for 
calculating  the  expected  return  from  complex  targets  certainly  has 
the  potential  to  be  a  useful  simulation  tool.  There  are  of  course 
questions  which  must  be  addressed  before  the  technique  could  be  used 
for  most  practical  applications.  Specifically,  can  practical  algorithms 
be  developed  which  are  assured  to  find  sdl  the  twinkles  and  initial 
specul2ir  points  for  a  complex  target?  What  happens  if  the  target 
surface  contains  discontinuities  or  is  made  up  of  sepzirate  bodies?  And 
what  about  the  effects  of  shadowing?  This  section  attempts  to  show 
that  each  of  these  questions  can  be  answered  in  a  positive  manner 
and  should  be  the  topics  of  further  research.  The  section  closes  with 
the  suggestion  that  such  further  research  edso  include  experimental 
validation  of  the  method  smd  comparison  with  other  suitable  methods. 

5.2  Efficiently  Searching  for  Twinkles  and  Initial  Specular  Points 

Developing  efficient  yet  thorough  edgorithms  for  finding  the  twinkles 
and  specular  points  associated  with  a  particular  trajectory  and  its 
beginning  is  an  important  next  step  in  the  development  of  the  overall 
simulation  method.  Once  the  twinkles  and  initial  specular  points  are 
located,  the  process  of  determining  the  motion  of  the  spec\ilar  points 
is  both  straightforward  and  efficient.  As  was  pointed  out  in  the  pre¬ 
vious  examples,  however,  if  either  a  twinkle  associated  with  a  birth  or 
an  initial  specular  point  is  overlooked,  it  is  possible  to  miss  an  entire 
specular  path.  Thus,  the  overall  simulation  method  will  be  no  more 
robust  than  the  method  for  finding  the  twinkles  and  initial  specular 
points.  This  section  begins  by  showing  two  fairly  obvious  methods 
for  reducing  the  area  of  the  target  which  must  be  searched  for  either 
twinkles  or  specular  points.  A  proposal  is  then  made  to  exploit  the 
variation-diminishing  property  of  B-spline  surfaces  to  develop  tech¬ 
niques  which  petition  individual  B-spline  patches  into  sections  which 
do  not  have  twinkles  or  specular  points  and  those  which  may  have. 

Regardless  of  the  actual  gradient  method  used  to  locally  search  for 
either  twinkles  or  specular  points,  it  would  be  profitable  to  develop 
techniques  which  quickly  eliminate  areas  to  be  searched:  i.e.,  areeis 
for  which  one  of  the  necessary  conditions  cannot  be  met.  When  the 
target  surface  is  described  by  a  mesh  of  B-spline  patches,  it  makes 
sense  to  consider  the  individual  patches  as  identities  to  be  tested. 


For  example,  patches  which  face  in  a  considerably  different  direction 
from  that  of  the  trajectory  origin  can  be  quickly  identified  and  need 
not  be  searched  at  all  when  initial  specular  points  are  sought.  The 
following  is  a  second  and  particularly  usefid  technique  for  eliminating 
large  numbers  of  patches  from  the  twinkle  search. 

Recall  that  if  the  function  l|G(w,  u,  A)||*  is  considered  to  be  the 
w- coordinate  in  an  orthogonal  u,v,w  coordinate  frame  as  was  shown 
in  figure  9,  the  condition  for  a  twinkle  implies  that  one  of  the  two 
principal  radii  of  curvature  of  the  newly  defined  surface  is  equid  to 
zero.  Let  vt^,  and  At„  be  the  vzdues  of  u,  v,  and  A  at  a  twin¬ 
kle.  Another  interpretation  is  that  there  exists  a  curve  on  the  surface 
passing  through  Vtu,  with  the  radius  of  curvature  originating  at 
the  trajectory  location,  a  fraction  At^  along  the  trajectory.  Clearly, 
if  a  particular  B-spline  patch  is  totally  convex  with  its  radius  point¬ 
ing  away  from  the  trajectory,  then  the  surface  (u,  v,  ||G(u,  v,  A)||*) 
will  also  be  convex  with  the  radius  pointing  away  from  the  trajec¬ 
tory.  And  therefore,  the  patch  cannot  contain  a  twinkle  (note  that 
this  does  not  imply  that  the  patch  will  not  have  specular  points).  In 
most  cases  a  particular  B-spline  patch  is  only  viewed  from  one  side 
(for  example  if  the  patch  is  part  of  a  closed  surface).  Therefore  such 
totally  convex  patches  can  be  eliminated  from  the  twinkle  search  for 
all  trajectories  of  interest.  Since  many  targets  consist  of  large  areas  of 
convex  cap  regions  (regions  where  the  curvature  does  not  change  sign 
and  points  toward  the  target  interior),  the  savings  from  the  reduced 
twinkle  search  can  be  substantisil.  This  technique  is  particularly  ef¬ 
ficient  since  the  patches  to  be  eliminated  are  found  only  once  for  a 
particular  target  and  their  identification  stored  away  for  rapid  future 
use  with  all  trajectories. 

Now  consider  the  problem  which  occurs  when  we  attempt  to  find 
sdl  the  specular  points  associated  with  the  initial  trajectory  location. 
After  those  patches  which  are  deemed  unable  to  contain  any  specular 
points  are  removed  from  the  search,  we  are  still  left  with  the  problem 
of  searching  individual  patches  for  an  unknown  number  of  specular 
points.  Since  specular  points  may  occur  arbitrarily  close  to  each  other, 
we  would  have  to  begin  specular  searches  arbitrarily  close  to  each 
other  to  ensure  that  all  the  points  are  found. 

Until  now,  the  second-order  continuity  of  the  B-spline  surface  is 
the  only  feature  of  the  surface  which  has  been  specifically  exploited. 
However,  an  important  property  of  B-spline  surfaces  is  their  variation- 
diminishing  property.  Loosely  stated,  the  variation-diminishing  prop¬ 
erty  says  that  each  bicubic  B-spline  surface  patch  on  a  uniform  rec- 


tilinear  mesh  will  have  no  more  undulations  than  the  16  three-di- 
mension2d  data  points  on  which  it  is  defined.  As  an  immediate  conse¬ 
quence,  if  we  wanted  to  find  local  maxima  (minima)  on  a  B-spline  tar¬ 
get,  it  would  only  be  necessary  to  search  areas  near  data  points  which 
are  themselves  local  maxima  (minima)  with  respect  to  neighboring 
data  points.  It  would  be  tempting  to  conclude  that  in  order  to  find 
all  the  locations  on  the  target  surface  for  which  the  distance  to  a  point 
on  the  trajectory  is  a  local  maximum  (minimum),  it  is  only  necessary 
to  search  near  data  points  for  which  the  distance  to  the  trajectory  is 
a  maximum  (minimum)  with  respect  to  its  neighboring  data  points. 
However,  such  a  conclusion  would  not  be  valid.  To  see  this,  it  is  only 
necessary  to  note  that  the  surface  (u,u,  ||G(u,v,  A)|p)  with  A  fixed 
is  a  higher  order  surface  than  the  original  B-sphne;  i.e.,  the  surface  is 
sixth  order  in  both  u  and  v  and  could  clearly  have  more  stationary 
points  than  the  original  bicubic  surface.  Fortunately,  the  resulting 
surface  is  still  just  the  tensor  product  of  two  sixth-order  polynomi¬ 
als,  and  may  be  represented  by  the  tensor  product  of  two  sixth-order 
B-splines.  The  resulting  higher  order  B-spline  surface  patch  will  be 
defined  by  a  set  of  49  points  (7x7  grid)  which  can  then  be  easily 
checked  for  possible  maxima  (minima).  The  effort  required  to  trans¬ 
form  the  distance  to  the  target  surface  into  the  higher-order  B-spline 
surface  may  be  considerable.  However,  ignoring  for  now  the  effects  of 
numerical  inaccuracies  in  the  computing  process,  it  should  be  possible, 
using  the  transformed  surface,  to  develop  ordinary  gradient  searches 
which  are  assured  to  converge  to  all  the  specular  points.  Since  the 
global  search  for  specular  points  need  only  be  conducted  once  for  each 
trajectory,  the  time  required  for  the  surface  transformation  may  be 
amply  rewarded  by  the  increase  in  robustness  of  the  overall  simulation 
method. 

It  may  of  course  be  possible  to  create  other  more  efficient  methods 
for  assuring  that  all  the  initial  specular  points  are  found.  At  the 
very  least,  the  method  described  above  can  serve  as  a  standard  for 
comparison. 

The  technique  of  using  the  variation-diminishing  property  of  B- 
spline  surfaces  to  ensure  that  all  the  specula  points  corresponding 
to  the  initial  trajectory  location  are  found  may  be  extended  to  give  a 
method  for  ensuring  that  all  the  twinkles  corresponding  to  a  particular 
trajectory  are  found.  Consider  the  following  function: 

M{u,v,A)  =  F^  +  Fi  +  F^, 

where  Fi,  F2,  and  F3  were  defined  in  equations  (2),  (3),  and  (7), 


39 


respectively.  Clearly,  M{u,v,  A)  is  greater  than  or  equal  to  zero  for  all 
u,  V,  and  w.  In  addition,  M{u*,v*,A*)  equal  to  zero  implies  that  u*, 
V*,  and  A*  are  the  parametric  coordinates  of  a  twinkle.  If  we  now  let 
M(u,v,  A),  u,  V,  and  A  be  the  orthogonal  axes  of  a  four-dimensional 
hyperspace,  then  twinkles  will  occur  for  values  of  u,  u,  and  A  such 
that  the  resulting  hypersurface  equals  zero  (which  also  happens  to 
correspond  to  a  local  minimum).  Analogous  to  the  previous  section, 
it  shoxild  be  possible  to  represent  the  resulting  hypersurface  by  a  three- 
way  tensor  product  of  appropriately  ordered  B-splines: 

Miu,v,A)  =  X:EE^o^i(u)B,(u)B,(A). 

i  j  k 

The  variation-diminishing  property  now  shows  that  twinkles  can  only 
occur  near  the  A,j/k  with  negative  values. 

Once  again,  this  observation  may  or  may  not  lead  to  the  best 
method  for  determining  the  number  and  locations  of  the  twinkles 
associated  with  a  particular  trajectory.  However,  it  does  indicate 
that  such  methods  do  exist  and  that  further  research  in  this  area  will 
probably  be  productive. 

5.3  Specular  Paths  as  Roots  of  Polynomial  Surfaces 

Consider  again  the  distance  surface  of  figure  9.  The  problem  of  find¬ 
ing  all  the  specular  points  associated  with  a  single  location  along  the 
trajectory  is  the  same  as  that  of  finding  all  the  extrema  (maxima  and 
minima)  of  the  surface.  If  the  surface  is  a  piecewise  bicubic  polyno¬ 
mial  surface,  as  it  would  be  if  the  original  target  were  a  piecewise 
bicubic  B-spline,  the  specular  points  are  just  the  simultaneous  roots 
of  the  surfaces  formed  by  taking  the  partied  derivatives  of  the  distance 
surface  with  respect  to  the  two  parameters  u  and  v.  In  this  context, 
twinkles  have  a  simple  intuitive  interpretation.  They  are  simply  the 
locations  at  which  the  imaginary  parts  of  the  complex  roots  go  from 
being  nonzero  to  zero,  or  vice  versa.  (If  this  is  not  immediately  obvi¬ 
ous,  just  consider  the  fact  that  any  tensor  product  polynomial  surface 
of  a  given  order  has  an  invariant  number  of  roots,  provided  we  allow 
the  roots  to  be  complex  and  count  multiple  roots.)  The  motion  of 
the  S/R  across  the  target  surface  results  in  a  continuous  distortion 
of  the  distance  surface.  As  such,  we  expect  that  roots  of  the  partial 
derivatives  of  the  surface  should  move  in  continuous  paths.  It  might 
appear  at  first  that  twinkles  represent  a  violation  of  the  continuity  of 
the  paths  of  the  roots.  However,  that  is  only  because  in  the  search 


for  specular  points,  we  were  only  interested  in  the  roots  with  zero 
imaginary  part.  Note  that  the  requirement  that  specular  points  must 
be  annihilated  or  created  in  pairs  can  be  seen  in  this  context  as  the 
fact  that  the  roots  which  have  nonzero  imaginary  parts  must  occur  in 
complex  conjugate  pairs. 

Consider  two  methods  for  tracking  specular  points  by  searching  for 
complex  veJues  of  the  u  and  v  (call  them  u  and  v)  parameters  which 
are  extrema  of  the  distance  surface.  The  first  is  a  brute-force  method. 
If  we  think  of  each  of  the  polynomial  B-spline  patches  as  being  defined 
over  the  entire  imaginary  plane,  then  each  patch  will  contain  the  same 
number  of  extrema  (of  course  some  of  the  extrema  may  be  repeated). 
It  should  be  possible  to  find  all  these  complex  extrema  for  a  single 
trajectory  location  and  for  each  patch,  then  track  the  motion  of  these 
complex  points  Jis  the  S/R  moves  along  the  trajectory.  Those  extrema 
which  have  zero  imaginary  parts  will  be  given  the  distinction  of  being 
specular  points.  Unfortunately,  if  the  surface  consists  of  thousands  of 
polynomial  patches,  then  even  more  thousands  of  complex  extrema 
will  have  to  be  tracked  for  all  positions  along  the  trajectory  to  assure 
the  invariance  in  the  number  of  extrema. 

For  a  second  method,  suppose  that  the  polynomial  patches  were 
linked  together  so  that  not  only  the  composite  surface  is  smooth  in 
the  sense  we  have  been  considering,  but  sJso  the  complex  extrema 
of  the  surface  are  continuous  across  patches.  With  such  a  surface, 
it  would  only  be  necessary  to  find  all  the  complex  extrema  with  real 
parts  in  the  range  for  which  we  consider  the  patch  to  describe  the 
surface;  i.e.,  the  real  parts  of  u  and  v  restricted  to  [0,1].  Then  as 
the  S/R  moves  along  the  trajectory,  the  only  way  that  an  extremum 
could  enter  a  new  patch  is  if  it  leaves  a  neighboring  patch  (ignoring 
for  the  moment  the  edges  of  the  composite  surface).  If  the  composite 
surface  is  closed,  finding  all  the  extrema  with  real  parts  in  [0,1]  will 
result  in  finding  an  invariant  set  of  complex  extrema  which  is  much 
reduced  from  that  of  the  previous  brute- force  method. 

This  method  has  the  potential  to  be  extremely  powerful  and  effi¬ 
cient.  Once  all  the  complex  extrema  are  found  for  any  location,  they 
may  be  tracked  to  the  beginning  of  any  desired  trajectory  and  then 
tracked  along  the  desired  trajectory.  Note  that  it  would  be  trivial  to 
include  arbitrary  nonlinear  trajectories.  Note  also  that  because  the 
invariant  number  of  extrema  need  only  be  determined  once  for  any 
given  target,  the  methods  used  to  find  the  extrema  can  be  extremely 
thorough. 

It  should  be  noted  that  the  search  for  complex  roots  C2ui  be  con- 
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siderably  more  difficult  than  the  search  for  strictly  real  roots.  One 
potentially  useful  method  for  finding  the  complex  roots  of  a  B-spline 
patch  would  be  to  find  all  the  roots  for  a  particulair  patch  (with  u  and 
V  extended  to  the  entire  imaginary  plane)  chosen  to  simplify  the  prob¬ 
lem.  Once  all  these  roots  are  found,  the  B-spline  coefficients  could  be 
modified  to  the  desired  B-spline  patch  while  the  complex  roots  are 
tracked.  This  method  would  assure  that  all  the  complex  roots  are 
found.  Since  neighboring  patches  tend  to  have  similar  coefficients, 
it  should  be  natural  and  efficient  to  continue  to  deform  patches  into 
neighboring  patches  while  tracking  the  effect  on  the  roots. 

Unfortunately,  cubic  B-spline  surfaces  as  used  in  this  report  do  not 
provide  continuity  of  the  complex  roots  across  patches.  One  would 
hope  that  just  as  raising  the  order  of  the  spline  allows  for  additional 
derivative  continuity,  perhaps  the  additional  freedom  of  raising  the 
spline  order  by  one  could  be  used  to  ensure  continuity  of  the  complex 
roots.  Because  of  the  very  powerful  potential  of  this  method,  the 
requirements  for  complex  root  continuity  of  the  B-spline  composite 
surface  should  be  investigated. 

Twinkle  Search  Through  Recursive  Subdivision 

It  was  pointed  out  in  section  5.2  that  if  a  portion  of  the  surface  is 
such  that  its  maximum  and  minimum  principal  radii  of  curvature  both 
point  toward  the  interior  of  the  surface,  then  no  twinkles  can  exist  on 
that  portion  of  the  surface  for  any  trajectory.  Consider  a  generaliza¬ 
tion  of  this  technique.  Each  point  on  a  patch  of  the  B-spline  quilt 
hats  associated  with  it  two  principal  radii  of  curvature  (the  maximum 
and  minimum  radii  of  curvature).  Suppose  the  maximum  and  min¬ 
imum  values  of  the  two  principal  radii  of  curvature  were  known  for 
the  entire  patch.  If  we  add  to  the  two  maxima  an  amount  to  account 
for  the  finite  extent  of  the  patch,  and  subtract  from  the  minima  the 
same  amount,  the  four  resulting  radii  will  provide  two  bounding  vol¬ 
umes,  one  of  which  must  be  intersected  by  the  trajectory  if  a  twinkle 
is  to  occur  on  that  patch;  i.e.,  each  principal  radius  of  curvature  will 
have  associated  with  it  a  range  guaranteed  to  include  its  maximum 
and  minimum  values  over  the  entire  patch.  If  the  trajectory  does  not 
pass  through  either  of  these  two  ranges,  then  the  patch  cannot  con¬ 
tain  a  twinkle.  Now  consider  what  happens  if  the  originzil  patch  is 
subdivided  into  four  subpatches.  The  maximum  and  minimum  val¬ 
ues  of  the  two  principal  radii  of  curvature  will  become  closer  to  the 
same  on  each  subpatch,  decreasing  the  volume  of  space  which  must 


be  intersected  by  the  trajectory  for  any  of  the  subpatches  to  contain 
a  twinkle  (the  extent  of  the  patch  wiU  decrease  as  well).  In  the  limit 
as  the  patch  is  recursively  subdivided,  two  concentric  spheres  will  be 
formed,  one  of  which  the  trajectory  must  intersect  in  order  for  the 
subpatch  to  contain  a  twinkle. 

The  above  method  can  be  turned  into  a  recursive  subdivision  tech¬ 
nique  for  finding  twinkles  by  adding  the  requirement  that  a  twinkle 
must  edso  be  a  specular  point.  Using  the  normal  at  the  center  of  a 
patch  along  with  the  maximum  curvature  for  the  patch,  we  can  cal¬ 
culate  a  “normal  cone”  which  includes  all  possible  surface  normals 
for  the  patch.  As  the  patch  is  subdivided,  each  cone  will  become 
smellier  in  extent  until,  in  the  limit,  the  cone  becomes  a  line.  Twin¬ 
kles  will  be  assumed  to  exist  only  on  those  subpatches  for  which  the 
trajectory,  normal  cone,  and  space  between  appropriate  concentric 
curvature  spheres  all  intersect.  The  subdivision  can  be  carried  out 
until  limited  by  the  precision  of  the  computer,  at  which  time  the 
remaining  candidate  subpatches  may  be  considered  the  locations  of 
twinkles.  Except  for  twinkles  that  lie  so  close  together  that  the  dif¬ 
ference  in  their  positions  czmnot  be  distinguished  with  the  precision 
of  the  computer,  the  recursive  subdivision  technique  should  result  in 
finding  all  the  twinkles  for  a  given  trajectory.  It  should  be  clear  to  the 
reader  that  the  normal  cone  described  above  could  also  be  used  by 
itself  to  define  a  recursive  subdivision  test  for  initial  specular  points. 

It  should  be  possible  to  simplify  the  search  for  the  maxima  and 
minima  of  the  two  principal  curvatures  on  each  subpatch  by  first  find¬ 
ing  all  stationary  points  of  the  principal  curvatures  for  each  original 
patch.  The  extrema  of  the  principal  curvatures  on  all  subpatches 
would  either  have  to  be  the  extrema  of  the  parent  patch  or  lie  along 
the  edge  of  the  subpatch,  greatly  simplifying  the  search  for  the  maxi¬ 
mum  and  minimum  curvatures  of  each  subpatch.  Note  that  the  search 
for  the  curvature  extrema  on  the  original  patches  would  only  have  to 
be  carried  out  once  for  the  model  since  curvature  extrema  are  indepen¬ 
dent  of  the  chosen  trajectory.  It  may  even  be  possible  to  find  efficient 
bounds  on  the  maximum  and  minimum  principal  curvatures  more  di¬ 
rectly  from  the  original  data  points  by  exploiting  the  properties  of  the 
B-spline  surface. 

5.5  Inclusion  of  Surface  Discontinuities 

The  surfaces  of  most  targets  of  interest  are  not  everywhere  second- 
order  continuous.  This  is  particularly  obvious  in  manmade  targets 
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which  generally  have  an  abundance  of  corners  and  edges.  While  it 
may  be  possible  to  model  edges  and  corners  as  second-order  continu¬ 
ous  surfaces  containing  locations  of  very  small  radius  of  curvature,  in 
many  czises  the  modeling  effort  will  be  made  much  easier  if  the  con¬ 
straint  of  second-order  continuity  everywhere  is  reduced;  that  is,  if 
the  surface  is  everywhere  second-order  continuous  except  along  edges 
or  corners. 

At  first  glance,  it  might  appear  that  the  use  of  twinkles  to  reduce 
the  computational  time  required  to  simulate  the  return  from  complex 
targets  wovild  be  of  little  use  when  the  surface  contziins  discontinuities 
of  either  its  first  or  second  derivatives;  in  such  a  situation,  it  would 
be  possible  for  specular  paths  to  begin  and  end  at  places  other  than 
twinkles.  Fortunately,  it  is  easy  to  see  that  for  all  the  target  surface 
which  is  second-order  continuous,  specular  paths  can  still  only  be  cre¬ 
ated  or  destroyed  at  twinkles.  Therefore,  once  the  twinkles  and  initial 
specular  points  are  found,  the  specular  paths  can  still  be  tracked,  as 
long  as  for  each  new  trajectory  location  all  the  discontinuities  are  also 
checked.  If  the  surface  discontinuities  occur  at  isolated  points  (cor¬ 
ners)  or  along  curves  (edges),  the  additional  search  required  by  the 
inclusion  of  discontinuities  will  be  reduced  in  order.  If  in  addition  the 
discontinuities  are  restricted  to  lie  on  parametric  curves  running  or¬ 
thogonal  to  either  of  the  two  surface  parameters,  the  problem  becomes 
particularly  simple.  For  example,  the  surface  normal  along  one  edge 
of  a  B-spline  patch  is  a  function  of  only  one  parameter  (with  proper 
choice  of  axes,  the  other  parameter  will  be  fixed  at  zero).  Thus  if  a 
discontinuity  lies  along  the  edge  of  a  B-spline  patch,  it  should  be  a 
straightforward  procedure  to  find  ^my  intersections  between  the  sur¬ 
face  normal  along  the  edge  and  the  trajectory.  The  intersections  could 
be  thought  of  as  “pseudo-twinkles”  since  they  represent  additional 
points  on  the  surface  at  which  a  specular  path  may  begin  or  end.  If  it 
is  necessary  to  include  discontinuities  along  parametric  curves  which 
are  not  orthogonal  to  either  of  the  two  surface  parameters,  the  algebra 
becomes  somewhat  more  unwieldy.  However,  it  should  still  be  a  fairly 
straightforward  task  to  find  the  pseudo-twinkles  along  the  curves  for 
an  entire  trajectory.  Thus  the  inclusion  of  curvature,  slope,  or  even 
position  discontinuities  in  the  target  surface  is  another  topic  whose 
investigation  should  prove  fruitful. 


5.6  Multiple  Bodies  and  Shadowing 

Backscatter  simulation  algorithms  based  on  twinkles  axe  just  as  ap¬ 
plicable  to  multibody  targets  as  they  are  to  single-body  targets.  Even 
for  multiple  bounces,  the  question  of  whether  the  bounces  occur  on 
separate  bodies  or  not  never  enters  into  the  problem.  The  problem 
with  finding  the  return  from  a  target  made  up  of  separate  bodies  is 
that  the  effect  of  shadowing  becomes  more  important;  i.e.,  it  is  more 
likely  that  the  ray  being  considered  will  have  passed  through  part  of 
the  target.  Recall  that  the  stationary  condition  for  multiple-bounce 
return  shown  earlier  (which  includes  the  single-bounce  ca.se)  does  not 
require  such  an  unobstructed  path. 

Questions  related  to  shadowing  are  almost  adways  handled  by  some 
form  of  ray  tracing.  However,  the  computational  savings  due  to  the 
use  of  twinkles  would  be  tot^llly  eliminated  if  it  were  still  necessjiry 
to  ray  trace  the  entire  scene  just  to  handle  the  shadowing.  Those 
specular  points  which  have  been  found  by  ignoring  potential  shadow¬ 
ing  we  refer  to  as  “potential  specular  points.”  Clearly,  the  problem 
of  determining  whether  or  not  a  potential  specular  point  is  an  actual 
specular  point  is  greatly  simplified  by  knowledge  of  its  position.  At 
most  it  will  be  necessary  to  check  only  the  paths  associated  with  each 
specular  point  for  intersections  with  the  target  surface  at  each  tra¬ 
jectory  location  to  determine  the  actual  “visibility”  of  the  return.  It 
may  even  be  possible  to  exploit  the  continuity  of  the  target  surface 
to  reduce  the  computational  time  reqmred  to  search  for  shadowing 
structures.  For  exsonple,  if  the  path  connecting  a  multiple-bounce 
return  is  to  go  from  being  unshadowed  to  shadowed,  some  portion 
of  the  path  must  become  osculatory  with  the  target  surface  at  some 
point.  And  of  course  the  same  must  hold  for  a  multiple  bounce  coming 
out  of  shadow.  Once  the  potential  specular  paths  have  been  deter¬ 
mined  (ignoring  potential  shadowing),  it  may  be  possible  to  exploit 
this  shadowing  constraint  to  rapidly  determine  which  portions  of  the 
paths  should  be  considered  actual  specular  paths  for  calctilating  the 
expected  return. 

5.7  Experimental  Validation  and  Comparison 

The  ultimate  goal  of  this  study  is  to  develop  a  useful  radar  backscat¬ 
ter  modeling  technique  for  such  targets  as  tanks,  aircraft,  and  ships. 
Clearly  then  an  experiment  of  interest  would  be  to  create  a  spline 
model  which  closely  resembles  the  surface  of  a  target  of  interest,  ex¬ 
perimentally  determine  the  backscatter  signatures  of  the  actual  target 
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for  many  trajectories,  and  compare  the  results  with  those  calculated 
from  the  spline  surface  using  a  geometric  optics  approach.  While 
such  experiments  will  eventually  determine  the  bottom-line  effective¬ 
ness  of  the  modeling  approach,  sm  intermediate  set  of  experiments 
using  a  target  developed  to  test  specific  assumptions  of  the  modeling 
method  should  also  prove  to  be  informative. 

Generally  speaking,  the  use  of  a  geometric  optics  assumption  for 
the  return  from  complex  targets  is  subject  to  two  broad  categories  of 
error.  The  first  type  of  error  is  due  to  the  surface  description  not  ad¬ 
equately  resembling  the  actual  target  surface  of  interest.  The  second 
category  of  error  is  due  to  the  geometric  optics  assumptions  them¬ 
selves  (i.e.,  everything  else).  For  a  target  such  as  a  tank,  there  are 
likely  to  be  a  number  of  items  which  do  not  fit  the  geometric  assump¬ 
tions  well.  Sharp  edges,  varying  surface  material  properties,  or  just 
plain  dirt  are  only  a  few  examples.  A  useful  experiment  therefore 
would  be  to  construct  a  complex  surface  which  satisfies  the  geometric 
optics  assumptions  to  a  high  degree  for  an  available  radar  frequency. 
The  surface  should  then  be  accurately  modeled  by  B-splines  and  the 
actual  and  simulated  returns  compared.  Comparisons  of  signatures 
using  various  frequencies  and  trajectories  should  lead  to  greater  in¬ 
sight  into  the  overall  modeling  problem  as  well  as  into  the  limits  of 
the  geometric  optics  approach. 

Once  a  geometric  optics  model  is  chosen  to  simulate  the  return 
from  complex  targets,  the  method  used  to  find  the  specular  points 
makes  no  diflference  in  the  calculation  of  the  return  (provided  each 
method  finds  all  the  specul2u:  points  for  each  location  of  the  trajec¬ 
tory).  However,  the  ability  to  rapidly  find  all  the  speculsir  points  as  a 
function  of  the  trajectory  location  does  give  the  ability  to  conduct  a 
more  extensive  validation  of  the  model.  Such  an  extensive  validation 
of  the  model  should  be  performed  for  a  complex  but  highly  controlled 
target  for  which  the  geometric  optics  assumptions  are  valid,  as  well 
as  for  actual  teirgets  of  interest. 
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6.  Summary 


Use  of  the  geometric  optics  approach  to  simulate  the  radar  backscatter 
from  composite  targets  cam  be  efficient  and  informative.  The  valid¬ 
ity  of  the  approach  depends  on  a  number  of  factors,  including  the 
target  surface  characteristics  and  the  wavelength  of  the  propagated 
signal.  The  use  of  B-splines  to  describe  complex  targets  for  use  in 
the  geometric  optics  model  is  an  obvious  extension  of  previous  meth¬ 
ods.  Unfortunately,  the  increased  complexity  of  the  model  means  am 
increaised  search  time  for  a  generally  unknown  number  of  speculair 
points.  Since  it  is  expected  that  the  speculair  points  must  be  found 
for  thousands  of  locations  along  each  trajectory,  ordinary  techniques 
which  search  the  entire  surface  for  specular  points,  at  each  location  on 
the  trajectory,  are  generally  too  computer  intensive  to  be  practical. 

Over  25  yeairs  ago,  Longuet-Higgins  [8]  pointed  out  the  continu¬ 
ity  of  specular  paths  on  smooth  surfaces  and  showed  the  conditions 
for  which  a  specular  path  may  become  discontinuous.  It  hais  been 
shown  in  this  paper  that  the  locations  of  specular  path  discontinu¬ 
ities  (known  as  twinkles)  may  be  predicted  for  cm  entire  trajectory. 
Once  the  locations  of  all  the  twinkles  associated  with  a  particular  tra¬ 
jectory  are  known,  it  is  only  necessary  to  track  the  motion  of  specular 
points  as  a  function  of  the  trajectory  location,  assured  that  creations 
and  annihilations  of  specular  paths  will  only  occur  at  the  twinkles.  It 
was  shown  that  the  method  extends  to  multiple  bounces  of  any  order. 

The  actual  algorithms  currently  used  to  find  twinkles  and  initial 
specular  points  are  implementations  of  simple  local  gradient  tech¬ 
niques,  which  offer  no  assurance  of  convergence  to  all  the  desired 
points.  The  shortcomings  of  these  techniques  were  discussed  and  sug¬ 
gestions  were  made  for  techniques  to  increase  the  robustness  of  the 
algorithms.  In  addition,  the  problems  of  surface  discontinuities  and 
shadowing  were  briefly  discussed,  with  the  conclusion  that  they  may 
be  included  in  the  model  without  negating  the  savings  incurred  by  us¬ 
ing  twinkles.  The  inclusion  of  surface  discontinuities  and  shadowing 
are  two  more  suggested  topics  of  further  research. 

Consider  the  problem  of  simulating  the  backscatter  from  a  target 
represented  by  the  equivalent  of  10,000  B-spline  surface  patches  for 
1000  locations  along  a  specified  trajectory.  The  author  anticipates 
that  pursuit  of  the  research  suggested  in  this  report  will  result  in 
algorithms  to  find  the  backscatter  signatures  from  such  encounters  in  a 
matter  of  minutes,  on  a  reasonably  high-powered  computer.  It  should 
be  possible  with  such  a  model  to  capture  the  significant  features  of 


most  complex  targets  of  interest.  Such  a  model  would  obviously  find 
important  applications  in  many  current  as  well  as  einticipated  projects 
using  high-frequency  radar  eis  a  source  of  information. 
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Appendix  A. 


Newton  Step  for  the  Location  of  a  Specular  Point 
on  a  B- Spline  Surface 

We  axe  interested  in  solving  the  matrix  equation 

/  \  r  Q  1  r  _  7?  ■ 

du  dv  1  -^1 

1  ^  ^  )  Si  ~  -Fi  ’ 


du  dv 


where 


a|lG(u,t,,A)|l^ 

i-i  :=  AJ  =  - ^ - , 


..  ai|G(u,u,A)||* 

If  we  let  the  x,y,z  coordinates  of  figure  2  in  the  body  of  the  report  be 
replaced  by  respectively,  the  equation  for  (|G(u,u,  A)j|^ 

becomes 

l|G(«,»,  A)||’  =  -  jcj)  +  x;  - 

»=1 

where  iZ,,  is  the  X*  component  of  the  surface  R(u,v).  Then, 

F.(u,.,A)  =  -2s|[A(x;-x;)  +  xi-fi,.(u,.)i5^^^J, 

=  -2E  {|a(x;  -  xj)  +  x;  - 
=  -2^  {lA(x;  -  x;)  +  x;  - 


-it  {|A(x;  -  x;)  +  x;  - 
»=1  \ 

r5iZ*<(u,u)  5i2*.(u,u)l  1 

^  ^  Jr 


=  -2t{lA(x;-x;)H 

r5iZ*<(u,u)  dRxT(u,v) 
[  du  dv 

dFi(u,v,A)  _  dFziUfVjA) 
dv  du 


_iEr_5£^ 


Appendix  A 


Let  Aij  be  the  vector  coefficients  ■  f  a  single  B-spline  patch  and  denote 
its  X,  y,  or  z  component  by  A*j-.  Then  for  the  1,  mth  patch  of  the 
bicubic  B-spline  surface  being  considered, 

<=i  i=i 

where  Bi[z)  for  i=l,  2,  3,  and  4  were  given  in  section  1.3  and  u,v  E 
[0,1].  Thus,  for  any  particular  patch  and  trajectory,  given  u,v,  and 
A,  the  Newton  step  [5i,52]  toward  a  specular  point  is  given  by 


If  u,v,  and  A  are  not  the  coordinates  of  a  twinkle,  then  the  inverse 
in  the  previous  equation  is  guaranteed  to  exist. 
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Appendix  B.  Newton  Step  for  the  Location  of  Twinkles 


We  are  interested  in  solving  the  matrix  equation 
/  dF^  dfi  ^  ^ 

du  dv  BA 

bF'j,  dF2  dF2 

Bn  dv  BA 

dF2  dF2  BF2 

^  du  dv  dA  ' 

where  F\  and  F^  are  the  same  as  in  appendix  A,  and 


-Fi] 

52 

= 

-F2 

53 

[-Fz 

Fz  :=  F3{u,v,A) 


dFi{u,v,A)  dF2{u,v,A) 
du  dv 


ldF,{u,v,A)y 


Expressions  for  Fi,  Fj,  and  ^  as  functions  of  u,v,  and 

A  were  found  in  appendix  A.  Expressions  for  the  remaining  terms  are 
as  follows: 


dFi{u,VyA)  ^ 

dA  ~  du  ’ 


dF2{u,v,A)  _  _  yi.dR,i{u,v) 

dA  ~  ^hx 

Dropping  the  u,v  and  A  arguments  for  brevity, 


du 

II 

d^F2\  (dF2\  fd^Fx] 

dudv  /  \  /  \  du^  ) 

-mi 

d^Fx\ 
dudv j  ’ 

dF3 

dv 

■(g)( 

d^F2\  (dF2\  (d^Fx\ 

dv^  j  \  dv  j  \5u9v  / 

-m( 

d^Fx\ 
dv^  }  ’ 

dFs 

dA 

■(£)( 

S‘F,\  (eF,\(S’F,\  /dFA 

dAdvj  \  y  \dAdu)  “  \  dv  J 

\  dAdv  J 

To  further  abbreviate  the  notation,  let  V’t 

=  [A(X‘  -  Jf 

I)  -xi- 

RA 

u,u)]: 

d^F2 

3 

,  a^F,.(u,v) 

dRgi(u,v) 

dudv 

-2  2- 

r  dudv^  ' 

dv  dudv 

du 

dv^  J 

a*Fi 

3 

0  V 

,  9®F,i(u,u) 

^dR,i(u,v)d^R^i(u,v) 

dRgi(u,v)  d*F*i(u,v)] 

dudv 

-2  2- 
»=1 

du^dv 

du  dudv 

du"^ 

dv  J 

d'^Fx 

3 

0  V 

1  d^R,*iu,v) 

ai?,.(u,v)  d*F,.(u,v) 

du^ 

-22- 

i=l 

du^ 

du  du^ 

) 

M 

— 

i 
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dv^ 

dAdv 

d^Fi 

dAdu 

d^Fi 

dAdv 


-2E 


i=l 

3 


V’i 


d^Rgi{u,v)  dRxi{u 


dv^ 


-3- 


dv 


1=1 


dv^ 


-2  2^(^2--Xi)  , 


1=1 

3 


i 

I 


,r)  d^R^i{u,v) 


-»  T  -  .1  •  ■  ^'  ■  a_.i'.‘  .  *U.»» 


K  '  i' 


Appendix  C.  Taylor  Coefficients  at  a  Twinkle  in  the  Rotated  Frame 

At  a  twinkle,  the  Taylor  coefficients  Oijk  are  given  by 


aijk  = 


du^dv^dA'^ 


Utw  ,Vtw 


where  Utto»  and  Atu,  are  the  coordinates  of  the  twinkle.  We  are 
interested  in  the  Taylor  coefficients  0300,  ooii,  Go2o>  ^2105  and  oioi  in 
the  u',  v',  A'  coordinate  system,  where 

u  =  u'  cos  $  —  V  sin  d  +  utm 

V  =  u'  sin  9  —  V  cos  9  +  Vty, . 

Let  Aj^,-  be  the  scalar  components  of  the  vector  coefficients  of  a 
single  B-spline  patch,  with  A**  representing  the  x,  y,  or  z  coordinates 
as  fc  =  1,  2,  or  3  respectively.  Then, 


l|G{u,»,A)||’  =  E 

fc=l 


<=i  j=i 


where  was  defined  in  appendix  B, 


^l|G| 

du' 


l|G||^ 

du'dA 

9^  l|G|l’ 

du'^ 


^|jG| 

a<i'» 


3 


=  -2E 


fc=i 

+  B, 

3 


ij=l 


4 

E 

».j=i 


Afc  » 


-2  E  -  Jff)  t  +  s.(“)4^ 


=  -2E 


*=i 


«J=1 


■^2 


-E 

».2=1  L 

3 

-2E 


^^i(«)A*o  /  A  ,  n  ruU* 


^2' 


k=l 


tj=l 


*  'S’BiW.* 


E 

t,i=i 


au'3 
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8°  l|G||’ 

du'^dv' 


5u'*  ‘•'  au  j 


3 

-3E 


£^a‘,b»  +  b<WA?>^ 


du"‘ 

aSiivV  * 


E 


M=i 


The  expressions  for  nrrd  ^  can  be  obtained  directly  from 

the  expressions  for  and  respectively  by  replacing  »'  with 

v'  throughout: 


-2E 


fc=l 


«ii=i 


4 

E^ 

•,;=i  L 


au'2  *'■  au' 
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a’fl^tr)  .s  „  r  t  ^  h‘ 

+  2  Ai^  2 

^n  r  tAt  8»B,(e)  8^^s  8^^  J  A  [SB.W  ,» 

+  B,(u)A,j-  +  Qy!  u 


-E 

<.j=i 


+a(»)At-^  E 

ij=l 


^®‘Wa?-B,(v)  +  2?^A‘?^ 


..i,d^Bj{v) 

+Bt(«)A‘-5;^ 


2E 


a.j-JV-'/  ■  -  Qyi  "‘j  ^u' 


. . , ,  au'  dv' 
,aB<(u),fcaB»]\ 


It  remains  to  determine  the  various  partial  derivatives  and 

First  note  that  |^  =  cos^,  0  =  -sin0,  ^  =  sin^,  and 

du  'dv  ^ 
dv 
dv* 


a  \a.»4  •  ***'''•  -  •  VV  ,  .  .  <*  J 

^  =  cose  and  further,  that  all  higher  order  derivatives  of  u  and 
Twith  respect  to  u'  or  »'  are  identically  zero.  It  is  easy  to  see  by 
induction  that 


a*Bi(u)  _  a*Bi(u) 
au* 


au'* 


du 


du' 


The  basis  step  is 


dBi(u)  dBi(u)  du  ,  dBi{u)  dv  _  dBjju)  du 

- =  r - -TT-.  T  5  Q.  /  ~  a..  a,,/ 


du' 


du  du'  dv  du'  du  du' 
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Appendix  C 


If  the  conjecture  holds  for  the  fcth  partial  derivative,  then 
d  d^B'iu) 


dv/  du'^ 


Similarly, 

d^Bijv)  _  ^Bi{v) 

dv'^  dv^ 

d^Biju)  _  d'‘Bi{u) 
dv'^  du^ 

d^Bijv)  _  d’^Bijv) 
du'^  dv'* 

And  further,  it  can  be  shown  that 


dv 


dv' 

du 


lk 


dv' 


dv  1 


du' 


du'^^dv'' 


a*+'5i(u) 

du 

k 

du 

du'*'^' 

du' 

dv' 

and 


dv 

k 

dv 

du>*dv’' 

dv'*'^' 

du' 

dv' 

d’^Biiu), 
du-  ‘ 

du 

du' 

fc-i 

d^u' 

du'\ 

d*+‘J3i(u) 

du 

du' 

du 

du' 

r  du 

■  k+l 

5u*+* 

1 

• 
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